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Abstract 

We consider the motion by mean curvature of an n-dimensional 
grapii over a time-dependent domain in M", intersecting K" at a con- 
stant angle. In the general case, we prove local existence for the cor- 
responding quasilinear parabolic equation with a free boundary, and 
derive a continuation criterion based on the second fundamental form. 
If the initial graph is concave, we show this is preserved, and that the 
solution exists only for finite time. This corresponds to a symmetric 
version of mean curvature motion of a network of hypersurfaces with 
triple junctions, with constant contact angle at the junctions. 



1. Time-dependent graphs with a contact angle condition. 

We consider a moving hypersurface in M"+^, with normal velocity 
equal to its mean curvature, assumed to be a graph over a time-dependent 
open set D{t) C M" (not necessarily bounded, or connected.) The (properly 
embedded) (n — l)-submanifold of intersection: 

T{t) = StnM" = dD{t) 

is a 'moving boundary'. Along T(t) we impose a constant-angle condition: 

(iV,e„+i)|r(f) = /3, 

where < /? < 1 is a constant and N is the upward unit normal of 'Ef 
'Mean curvature motion' is defined by the law: 

Vn = H, 

where V/v = {V, N), with V = dfF the velocity vector in a given parametriza- 
tion F{t) of Tit {V depends on the parametrization, while V/v does not). A 
particular parametrization yields 'mean curvature flow': 

dtF = HN. 

For graphs, it is natural to consider 'graph mean curvature motion': if = 
graph w{t) for a function w{t) : D{t) — > M, imposing {dtF,N) = H with 
F{y,t) = [y,w{y,t)] for y G D{t), we find: 



wt = + \Dw\'^H 
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(and the velocity is vertical, dtF = wten+i)- With the contact angle condi- 
tion, we obtain a free boundary problem for a quasilinear PDE: 

( wt = g'^{Dw )wy in D {t), 

\ w = 0, P^l + \Dw\'^ = 1 on dD{t), 

where g^^{Dw) = — WiWj/{l + is the inverse metric matrix. 

Remark 1.1. It is easy to see that the constant-angle boundary condition 
is incompatible with mean curvature flow parametrized over a fixed domain 
Dq: on ODq we would have {F,en+i) = 0, leading to {dtF,en+i) = 0, 
incompatible with dtF = HN and (A'^, e„+i) = If we parametrize over a 
time-dependent domain, mean curvature flow leads to a normal velocity for 
the moving boundary that is difficult to control; hence we chose to analyze 
the geometry of the motion in terms of the 'graph m.c.m.' parametrization. 

To establish short-time existence (in parabolic Holder spaces) we will 
work with a third parametrization of the motion, defined over a fixed do- 
main: 

F{t) : Do M"+\ F{x, t) = [ip{x, t), u{x, t)] G x M, 

where (p{t) : Dq D{t) is a diffeomorphism and F is a solution of the 
parabolic system: 

Ft = g'^{DF)F,j, 

where gij = {Fi,Fj) is the induced metric on and g^^ is the inverse metric 
matrix. 

In the first part of the paper (sections 3 to 8) we prove the following 
short-time existence theorem (on Q := Dq x [0, T]): 

Theorem 1.1. Let Sq C M"+^ be a C^+" graph over Dq C M" satisfying 
the contact and angle conditions at ODq. There exists a parametrization 

-^0 = ['/'Oj'Wo] £ C'^"'""(Z?o) of Eq (where a = a"^ and cpo is a diffeomorphism 
of Dq satisfying the 'orthogonality conditions' on ODq), T > depending 
only on Sq and a unique solution F G C^^""'-^"'""/^(Q-^; of the system: 

r dtF - g'i {DF)didjF = 0, F G M" x R 

\ =0, N^+\DF)\Qj,^=P, 

with initial data Fq and where ip{t) : Dq D{t) C is a diffeomorphism 
satisfying, in addition, the 'orthogonality conditions' at ODq (described in 
section 3.) 
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The system and boundary conditions are discussed in more detail in sec- 
tion 3. Sections 4, 5, and 6 deal with compatibility at t = 0, linearization 
and the verification that the boundary conditions satisfy 'complementarity'. 
In particular, adjusting the initial diffeomorphism ipQ to ensure compatibil- 
ity (section 4) leads to the 'loss of differentiability' seen in theorem 1.1. The 
required estimates in Holder spaces for the linearized system are described 
in section 7, and the proof concluded (by a fixed-point argument) in section 
8. While the general scheme is standard, since we are dealing with a free 
boundary problem with somewhat non-standard boundary conditions, de- 
tails are included. Free boundary-type problems for mean curvature motion 
of graphs have apparently not previously been considered. 

We describe the evolution equations in the rotationally symmetric case 
in section 9 (including a stationary example for the exterior problem) and 
the extension to the case of a graph motion Sj intersecting fixed support 
hypersurfaces orthogonally in section 10. 

The original motivation for this work was to establish (by classical parabolic 
PDE methods) existence-uniqueness for mean curvature motion of networks 
of surfaces meeting along triple junctions with constant-angle conditions. 
One can use a motion of graphs with constant contact angle to produce 
examples of 'triple junction motion': three hypersurfaces moving by mean 
curvature meeting along an (n — l)-dimensional submanifold $](t) so that 
the three normals make constant angles (say, 120 degrees) along T{t). We 
simply reflect on M", so the hypersurfaces are and M" — D{t). If 

S( = graph w{t) with > 0, the system is embedded in R""*"^. This is mean 
curvature motion of a 'symmetric triple junction of graphs'. 

Short-time existence holds for general triple junctions of graphs moving 
by mean curvature with constant 120 degree angles at the junction, pro- 
vided a compatibility condition holds along the junction (see section 16). 
Since the free-boundary problem is easier to understand in the symmetric 
case, we decided to do this first. In addition, in the present case it is pos- 
sible to go further towards a geometric global existence result. Motivated 
by recent work on 'lens-type' curve networks [7], in the second part of the 
paper (sections 11-15) we consider continuation criteria and preservation of 
concavity. Since we chose to develop these results for graph motion with a 
free boundary, although the general lines of proof (via maximum principles) 
have precedents, the details of the arguments are new. For example, section 
12 contains an extension of the maximum principle for symmetric tensors 
with Neumann- type boundary conditions given in 113], which in our setting 
allows one to show preservation of weak concavity in general. Section 14 
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includes a continuation criterion for the flow. The results obtained in sec- 
tions 11-15 are summarized in the following theorem, (h denotes the second 
fundamental form, pulled back to a symmetric 2-tensor on D{t).) 

Theorem 1.2. If Sq is weakly concave (/i < at t = 0), this is preserved 
by the evolution. Let T^ax be the maximal existence time for the evolution. 
If the mean curvature of Eq is strictly negative (supsg H = Hq < 0), then 
Tmax is finite. Assuming Tmax < co, we have: 

limsup[sup(|/i|g + |V*""/i*'^"|g)] = oo 

(if n = 2, in the concave case). If there is no 'gradient blowup' at Tmax, the 
hypersurface contracts to a compact convex subset of M" as t — > T^ax ■ 

Remark 1.2. We have not yet proved that the diameter tends to zero 
as t ^ Tmax, though this seems likely based on the experience with curves 
[7] (in the absence of 'gradient blowup'.) It is an interesting question (even 
in the concave case, for n = 2) whether the latter can really occur (that is, 
suppj I V*""/i*'*"|g ^ oo as t — > Tmax, while \h\g remains bounded on F^.) 

Acknowledgments. It is a pleasure to thank Nicholas Alikakos for origi- 
nally proposing to consider the problem of mean curvature flow for networks 
of surfaces meeting at constant angles, and for his interest in this work. Most 
of the work on short-time existence was done during a stay at the Max- 
Planck Institute for Gravitational Physics in Golm (January- June, 2007); I 
am grateful to the Max-Planck Society for supporting the visit and to the 
Geometric Analysis group at the AEI for the invitation. Finally, thanks to 
Mariel Saez for communicating the results of the Lens Seminar ( [7J ) and of 
her recent work on mean curvature flow of networks (partly in collaboration 
with Rafe Mazzeo, |11]). 

2. Normal velocity of the moving boundary. The evolution is nat- 
urally supplied with initial data Sq, a graph meeting R""*"^ at the prescribed 
angle. Since we are interested in classical solutions in the parabolic Holder 
space we expect an additional compatibility condition at t = 0. 

We discuss this first for graph m.c.m. w{y,t). 

Denote by T(t) a global parametrization of dD{t) (with domain in a 
fixed manifold, and 'space variables' left implicit). Differentiating in t the 
'contact condition' w{T{t),t) = 0, we find: 

wt + {Dw,t{t)) =0. 



4 



Denote by n = the unit normal vector field to r(t), chosen so that the 
directional derivative dnW > 0. The contact condition also implies the gra- 
dient of w is purely normal: 

Dw\QD(t) = {dnw)n. 

Combining this with the angle condition, and bearing in mind that dnW\Y(t) > 
0, we find: 

dnW =jOn dD{t), l3o := V^'P^- 

(In fact, this is a more convenient form of the 'angle' boundary condition 
for w, since it is linear.) Thus, on dD{t): 

= Vl + {dnWfH = Wt = -{t{t),n)dnW = -t'n{t)j, 

and we find the normal velocity of the moving boundary (independent of 
the parametrization of F^): 

rn = -^%(t), 

which in particular must hold at t = 0. Note that we don't get a 'com- 
patibility condition' in the usual sense (of a constraint on the 2-jet of the 
initial data), but instead an equation of motion for the moving boundary. 
(Later, in the fixed-domain formulation, we will have to deal with a real 
compatibility condition) . 

Remark 2.1. We remark that for more general (non-symmetric, non-flat) 
triple junctions with 120 degree angles, the condition: 

H^+H^ = ij3 on r(^^ 

must hold at the junction (for graphs, oriented by the upward normal), which 
in particular gives a geometric constraint on the initial data, for classical 
evolution in C^+^'i+^Z^. This is automatic in the symmetric case {w = 
—w^,w^ = 0), since = and = trgid'^w^ for / = 1, 2. 

3. Choice of 'gauge'. It is traditional in moving boundary problems 
to parametrize the time-dependent domain D(t) of the unknown w{y,t) by 
a time-dependent diffeomorphism: 

y = ^{x,t), ip{t):Do^D{t), 
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and then derive the equation satisfied by the coordinate-changed function 
from the equation for w (see e.g. |9] or {14j). Motivated by the work on 
curve networks (|10j) we will, instead, consider a general parametrization: 

F : DqX [0, T] M"+\ F{x, t) = [ip{x, t),u{x, t)] G x M 

and derive an equation for F directly from the definition of mean curvature 
motion: 

{dtF,N) = H. 

(We'll still assume ip{t) : Dq — > D{t) is a diffeomorphism.) The first and 
second fundamental forms are given by: 

gij = {Fi ,Fj), AiFi ,Fj) = {Fi^ , N) . 

(Notation: DF = Fid, D'^F{ei, Cj) = Fij, (a) is the standard basis of M"+^) 
The mean curvature is the trace of A in the induced metric: 

H = {g'\DF)Fi,,N). 

The equation for F is: 

{dtF-g'^{DF)Fij,N) =0. 

There is a natural 'gauge choice' yielding a quasilinear parabolic system: 

dtF-g'^{DF)Fij=0. 

We will sometimes refer to this as the 'split gauge', since in terms of the 
components F = [ip, u] we have the essentially decoupled system: 

f dtu- g'^{Dif,Du)uij = 0, 
\ dtip - g'^{Dip,Du)ipij = 0. 

The splitting is useful to state the boundary conditions: 

/ ''^\dDo—^ ('contact'), 

\ N^+\D^,Du)\aD,=P ('angle'). 

We immediately see there is a problem, since we have 2 scalar boundary 
conditions for n + 1 unknowns (and no moving boundary to help!) Our 
solution to this is to introduce n — 1 additional 'orthogonality conditions' at 
the boundary for the parametrization '^{t). We impose: 

{Drip, Dn(p)\dDo = ('orthogonality') , 
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for any r G TODq, where n denotes the inward unit normal to Dq. (We fix 
a tubular neighborhood J\f of ODq and extend n to so that dnii = in 
Af.) 

Geometrically, the 'orthogonality' boundary condition has precedent in 
a method often adopted when dealing with the evolution of hypersurfaces 
in M"+^ intersecting a fixed n-dimensional 'support surface' orthogonally 
(see e.g. [15]): one replaces vanishing inner product of the unit normals 
(a single scalar condition) by a stronger Neumann-type condition for the 
parametrization, corresponding to n — 1 scalar conditions. (More details are 
given in Section 10.) 

The system must also be supplied with initial data. We assume given an 
initial hypersurface Sq, the graph of a C^"^" function Uo{x) defined in the 
(j3+a (domain Dq C M". (The reason for this choice of differentiability class 
will be seen later.) It would seem natural to set (po = Idoo, ^^^^ causes 
problems (related to compatibility; see Section 4 below). We do require the 
1-jet of ipo at the boundary to be that of the identity: 

^0\dDo = ^d, DipoiQo = I. 

(In particular, the orthogonality condition holds at t = 0.) 

We need a more explicit expression for the unit normal, and for that we 
use the 'vector product': 



N{D^,Du) := (-l)"det 



ei • • • Cn+l 
DF^ ■ ■ ■ DF""-^^ 



Dif^ ... L>(^" Du 



:= [J{Dip,Du),J^] G M" X M, 

where DF^ G for i = 1, . . . n + 1, J,^ > is the jacobian of and (— 1)" 
is introduced to make sure the last component is positive. J{Dip, Du) is an 
M^-valued multilinear form, linear in the components Uj of Du and of weight 
n — 1 in the components of Dip. It is easy to check that J(I, Du) = —Du. 
The unit normal is: 

N{Dip,Du) = N{Dip,Du)/{\JiD^,Du)\'^ + (J^)^)^/^. 

Thus the angle condition may be stated in the form: 

l3[\J{Du,Dip)f + (J^)2]f/Jo = "^^I^Do' 
and we lose nothing by squaring it: 

B{Dip,Du) := (3^\J{Du,D^)\' - P^Wfen, = 0. 
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4. Compatibility and the choice of ipQ. Assume -D(/3o|a£)o ~ ^' 
Differentiating in t tlie contact condition u\q£,q = and evaluating at t = 0, 
we find: 

= g"^ {I, DuQ)uQij = gl^uoij on ODq. 
To interpret tliis condition, consider the mean curvature at t = 0, on dDo: 

Ho = —[{J{I,Duo),gQ(poij) + Jipo9oUoij], 

where: 

vo = [\J{I,Duo)\' + J^jf/J^ = {\Duo\' + = ^, 

using (recall /3o := \/l — /3^): 

J(I,Duo) = -Duq = -{DnUo)n = -^n 
on ODq. Thus the compatibility condition is equivalent to: 



Ho\dDo = -Po9o{^Oij,n)\ 



dDo- 



This implies we can't choose (/?o = Id (on all of -Do); unless ^fo|apo = 0, a 
constraint not present in the geometric problem (as seen above) |^ Instead, 
regarding Hq as given (by Sq), and using: 



we find the compatibility constraint: 



{{6ij - (3on'n^)(poij,n) = -^Ho on ODq. 
Given the zero and first order constraints on (po, this can also be written as: 

n''n^(poij,n) = --J-^Hq on ODq. 

P PO 

The next lemma shows this can be solved. 

^The compatibility condition Ho\gDa ~ does occur for graph m.c.m. with Dirichlet 
boundary conditions in a mean-convex domain |6j. 
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Lemma 4.1. Let Dq C M" be a uniformly C^"*""^ domain (possibly 
unbounded), h G C°(5L>o) (0 < a < 1). 

(i) One may find a diffeomorphism G Dif ^^"'{Dq) satisfying on ^Dq- 

ip = Id, d(f = 1, n ■ cP'Lp{n, n) = h. 

(ii) More generally, given a non-vanishing vector field e G C"^+°((9Doi K"), 
one may find G Dif /'^^"'{Dq) satisfying on SZ^q- 

(f = Id, dnf = e, n ■ d^(p{n, n) = h. 

If ODq has two components, we may even require ip to satisfy the conditions 
in parts (i) and (ii) at OiDq, 820^ (resp.), with different functions h. (This 
will be needed in section 10). 

As usual, a domain is 'uniformly C^+"' if at each boundary point there 
are local charts to the upper half-space (of class C^^°), defined on balls of 
uniform radius, and with uniform bounds on the C^'*'" norms of the charts 
and their inverses. 

Remarks: 4.1. The proof is given in Appendix 1. 

4.2. Note that, in particular, ip satisfies the orthogonality conditions at 

4.3. It is at this step in the proof that we have a drop in regularity: 
for C'^^" local solutions, we require C'^"^" initial data. While this is not 
unexpected in free-boundary problems (see e.g. [9j), I don't know a coun- 
terexample to the lemma if Dq is assumed to be a C^^° domain. 

4.4. In our application of the lemma, we in fact have h G C^'^'^(9-Do)) 
but this does not imply higher regularity for (p. 

5. Linearization. The evolution equation and boundary conditions in 
'split gauge' are: 

' Ft-g'^{DF)Fij = 0, 

U\dDo = 0; 

B{D^,Du)\QDo = 0' 

where: 

0{D^) := {D^ip,Dnip). 

Here D^ip = Dtp — {dn^){-,n) is an M"-valued (n — l)-form on ODq. We'll 

prove short-time existence for this system (with initial data 

by the usual fixed-point argument based on linear parabolic theory. Given 
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F = [(p,u\ in a suitable ball in this Holder space with center Fq = [(pq^uq], 
it suffices to consider the 'pseudolinearization' of the system: 

Ft - g'^{DFo)Fij = [g''{DF) - g'^{DFo)]Fj := T{F, Fq) := T (LPDE); 

a fixed point of the map F ^ F corresponds to a solution of the quasilinear 
equation. 

For the nonlinear boundary conditions, we need the honest lineariza- 
tion at Fq. For the angle condition, a computation using the boundary 
constraints on uq and ipQ yields: 

^CoB[Dip, Du] = -PPodnU - Pi{dn(p, n). 
The corresponding linear boundary condition will be: 

(3PodnU + Piidnif, n) = B{DF, DFq) := B, 

where: 

2B{DF\DF^) := B{Dipi, Dui)-B{Dip2, Du2)-CoB[D{^pi-ip2), D{ui-U2)], 
and we used: 

~jCQ[Dipo,Duo]\dDo = PPodnUQ + Po{dn'PO,n)\dDo = 0- 

Also, B{D{pQ, DuQ)^g£)^^ = 0, so at a fixed point B{Dip,Du)\QD^^ = 0. 

Linearizing the orthogonality boundary condition, we find that CQO[Dip\ 
is the (n — l)-form on dD^: 

(with sum over repeated indices.) The corresponding linear boundary con- 
dition is: 

{dn^p^prof) + {D'^Lp,n) = -Q{Dp,D(po) ■= ^, 
where proj^ denotes orthogonal projection M" — > TODq, and: 

n{D^i,Dip2) := 0{Dipi) - 0{Dip2) - CoOlDipi - Dip2], 

and we used: 

CoO[Dipo]\Qjj^ = {{dnipofr) + {D'^^o,n)\aDo = 0- 
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6. Complementarity. We wish to apply linear existence theory to the 
system: 

Ft-g'^{DFo)Fij=P, 
with boundary conditions at dD^: 



{LBC) { 
and initial conditions: 



(3(3odnU + (3q {dn(p, n) = 8, 



ut=o = uo, ipt=o = ipo- 

It is easy to see that the initial data satisfy the linearized boundary condi- 
tions, and above we constructed (po so as to guarantee g^^ {Duq, DipQ)iiQij\^Qj)^ = 
0. (There is no first-order compatibility condition for ipQ.) Thus the linear 
system satisfies the required compatibility at t = 0. 

Since the linearized boundary conditions are slightly non-standard, we 
must verify they satisfy the 'complementarity' (Lopatinski-Shapiro) con- 
ditions. We fix xq G ODq and introduce adapted coordinates (p, <j) in a 
neighborhood A/q C of xq in Dq: 

X eJ\fo^ X = To{a) + pn{a), a = {aa) € U, 

where Tq : W ^ M" is a local chart for ODq at xq {U C R""^ open). This 
defines a basis of tangential vector fields in ro(^), and we may assume that, 
at Xq: {Ta,Tb) = Sab and VraTb{xo) = (for the induced connection on 
TdDo). Let U and tj; be defined in (-pi,0) xU x [0,T] by: 

U{p, (7, t) = u{To (a) + pn{a) , t) , ^(p, a, t) = ip{To (a) + pn{a) , t) . 

In these coordinates, the induced metric is written (in 'block form'): 



[9{DFo)] = 



J |t=0 







In- 



at t = and xq. 
We have: 



Upp = D^u{n,n) (since dnfi = 0), 



11 



Uab = D^u{Ta, n) + Du ■ V TaU = D^u{Ta, u) at Xq, 

and we don't need Upa, since gpa = at xq. 
Thus: 

trg,D''u{xQ) = P^D\{n,n)+J2DMra,ra) = 0'Upp+Y,Uaa ■= P'^Upp+A^U, 

a a 

and, likewise: 

trg^D^ifixo) = 0^il)pp + A^V- 
For the hnearized orthogonahty operator, note that, at xq: 

CqO[DiI,] = i{^p,Ta)+i^a,n))Ta. 

Putting everything together, the hnear system to consider at xo is: 



Ut - f3^Upp - KU = 0, 



with boundary conditions: 

U\p=o = 0, 

Poi^'pjn) + PUp\p=o = b{a,t), 

{i^p,ra) + {^Pa,n)\p=o = UJa{(T,t), a = l,...n- 1. 

Now take Fourier transform in cr G M"~^ (corresponding to ^ G M"~^), 
Laplace transform in t (corresponding to p G C) to obtain: 

U{p,^,p) G C, Vi(p,C,p) G ^ G W\p G C,p < 0. 

In transformed variables, we obtain the system of linear ODE (in p < 0, for 
fixed (^,p)): 

li''Upp-{p+mU = Q, 
Writing the solution in the form: 







- C7(0) " 









we find the characteristic equation /3^7^ + p+ = 0, and choose the root 
7 so that i7 = (1//3)\/A (where A = p + and we take the branch of 
Re{y/K) > 0). Here (p,^ G A, where: 

A = {{p,0 eCx R"-i; \p\ + lel > 0,Re{p) > 
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Thus the solutions decay as p — > — oo. Let W+ be the space of such decaying 
solutions, dimcW^ = n — 1. The relevant boundary operator on is: 

U{0) 

(a vector in C x C x C""^). 

The 'complementarity condition' (see e.g. [3]) is the statement that 
B is a linear isomorphism from to C""^^. With respect to the basis 
{^7(0), (Vi(0),n), {ip{0),Ta)} of W+, the matrix of B is (in 'block form'): 

1 [0]ix(„-i) " 

-^/A -f^/A [0]ix(„-i) . 

[0](n-l)xl Ka](n- l)xl "^In-l _ 

This is triangular with non-zero diagonal entries for every (p, ^) G A. Hence 
B is an isomorphism. 

7. Estimates in Holder spaces. 

For the fixed-point argument based on the linear system, we need esti- 
mates for ||.F||q,, I of two types: 'mapping' and 'contraction' 
estimates. 

A bit more precisely, for T > 0, R > and Q-^ = Dq x [0, T] consider 
the open ball: 

Bl = {Fe C2+"'i+"/2(Q^,M"+i); ||F - Folb+a < R,F\t=o = Fq}. 

{Fq = [ipo, uq] is defined from the initial surface Sq, via Lemma 4.1.) Solving 
the linear system with 'right-hand side' defined by -F G BJ^ defines a map 
¥ : F F, and we need to verify that, for suitable choices of T and R, ¥ 
maps into BJ^ and is a contraction. 

Remark: The argument that follows is standard, and the experienced 
reader may want to skip to the statement of local existence at the end 
of the next section. On the other hand, the result is not covered by any 
general theorem proved in detail in a reference known to the author, and 
some readers may find it useful to have all the details included. Another 
reason is that, although the 'right hand sides' are clearly quadratic, without 
explicit expressions one might run into trouble with compositions (which 



U 



U 



I „ — n 
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behave poorly in Holder spaces), or when appealing to 'Taylor remainder 
arguments' if the domain is not convex. 

For 'mapping', we need estimates of the form: 

\\J^{F,Fo)\\a + \\BiDF,DFo)\\i+a + \\n{D^,Dipo)\\i+o. decays as T ^ 0+, 
and for 'contraction': 

mF\F'')-J^{F^,F')\\^MmDF\DF^)\U+^+\\n{D<p\D 

where fi{T) ^ as T ^ 0+. 

Notation: The {a, a/2) norms arc taken on , the (1 + a, (1 + a)/2) 
norms on dDo x [0, T]). Double bars without an index refer to the (2 + a, 1 + 
a/2) norm, single bars to supremum norms over Q^, and parabolic norms 
are indexed by their spatial regularity (q for (a, a/2), etc.) In general, we 
use brackets for Holder-type difference quotients. 

We deal with the estimates for the 'forcing term' first. Consider the 
map 

g : Imm{W,W+'^) GLn 

which associates to the linear immersion A the inverse matrix of {{Ai, Aj))f^-^, 
inner products of the rows of A. Q is smooth, in particular locally Lipschitz 
in the space W of linear immersions. Hence if F^,F'^ are maps — > M""*"^, 
such that DF' G C°^'°/2(Q'r) and DF'{z) G K for all z e , where K CW 
is a fixed compact set, we have the bound: 

\\g{DF') - g{DF')\\a < ck\\D{F' - f2)||„. 

In fact our maps are in (72+a,i+Q!/2^ jjpi ^ (ji+a,^ ^ From this higher 
regularity we obtain the decay as T — > 0+. Assuming F^\t=o = F'^\t=o, we 
have: 

\D{F^ - F2)| < [D{F^ - F^)]f^h^. 

Now recall the elementary fact: if D C is a uniformly domain 
(not necessarily convex or bounded), and / G C^{D) with a G (0,1), we 
have: [/](") < Cd||/||ci. (Here 'uniformly C^' means D can be covered by 
countably many balls of a fixed radius, which are domains of manifold- 
with-boundary local charts for D, with uniform bounds for the charts 
and their inverses. The constant Cd depends on those bounds.) Applying 
this to DF, where F = F^ — F"^ vanishes identically at f = 0, and assuming 
r < 1: 
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[DF\^^^ <c{\DF\+\D'^F\)<c{[DF]\ ' ^T^-^ + [D'^F]\^'t"/'^)<c\\F\\T'^/^, 
(where c depends on Dq) and similarly for the oscillation in t: 
[DF\[^'^ < [L>F]J^Vi/2 < ||F||T^/^ 

so we have: 

\\D{F^ - F^)\\^ < c\\F'^ - F^llr"/^. 
We conclude, under the assumption F^ = F^ at t = 0: 

\\g{DF^) - g{DF^)\\a < ck\\F^ - F^\\T^/\ 

In particular, applying this to F and Fq, we find: 

WiGiDF) - g{DFo))D^F\\a < ck\\F - Fo||r"/2||F||, 

and for F^ and F^ coinciding at t = 0: 

\\{g{DF^) - g{DF^))D^F^\\a < ck\\F^ - F'^\\T''/^\\F^\\, 

as well as: 

\\{g{DF^) - g{DFo)){D'^F^ - D'^F'^)\\a < ck\\F^ - Fo\\T''/'^\\F^ - F\ 
so we have the mapping and contraction estimates for J^{F, Fq) and J^{F^ , Fq) — 

Lemma 7.1. Assume F,Fq,F^,F'^ are in C2+°'i+"/2(Q^; R^+i) and 
have the same initial values, and that DF, DFq, DF^ , DF"^ all take values 
in the compact subset K of Imm{W^ ,W^~^^) . Then: 

||^(^,i^0)||a<CK||^-i^0||||^||r"/2, 

mF\Fo) - J^{F^,Fo)\\a < ck{\\F'\\ + \\F^- Fo\\)T^/^\\F^ - F\ 
In particular, if F 6 BJ^: 

mF,Fo)\\a<coRT''/^. 

If F^,F^ e Sg, we have: 

mF\Fo)-HF^Fo)\\a < coT"/2||Fi - F% 
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(The constant cq depends only on the data at t = 0, and we assume T < 1, 
R < 1). 

Turning to the orthogonahty boundary condition, first observe that: 

which has quadratic structure. Using a local frame (Ta)^zJ for TODq, we 
find the components ^la- 

(summation convention, i, j = 1, . . . n), so is a sum of terms of the form: 
b{x)D{ip^ — ip'^)D{ip^ — ip*), where b{x) = n^rl and the (p^ coincide at i = 0. 
It is then not hard to show that: 

\\b{x)D{<p' - ^2)D(^3 _ ^4)||^_^^ < C||6||i+„||^1 - ^2||||^3 _ ^4||2.a 

with c depending on the norms of local charts for Dq. To bound the norm 

||n<8)ra||i+a, note \n\\Ta\ < 1, \D{n^Ta)\ < \Dn\ + \DTa\ and [D(n0ra)]l"^ < 
[Dn]!"-' + [DtoIx- Since n = —{P/Po)Duo on ODq (and ODq is a level set of 
uq), we clearly have: 

\\Dn\\a + ll-DTalla < c||i:)^Mo||a < cH^UoH- 

We summarize the conclusion in the following lemma. 

Lemma 7.2. Assume (p, ipo, (p^j^p"^ G C^+"'^+"/^((5^; M") have the same 
initial values. Then: 

\2rpa 



||0(D(^,D<^o)||i+a < co||no||||(^-(^o|rr'^ 
and 

\\niDip\D^^)\\i+^ < collt/olKllv'' - M\ + W^^ - ^o||)T"||(^i - ip\ 

with Co depending only on the data at t = 0. In particular, if F = [(p, u] € 
B-^, we have: 

\\n{D<p,D^o)\\i+a<coR''T^, 
and for = [(p^ , u^] e D^, I = 1, 2: 

\\n{D<p\D^^)\\i+a < coRT"\\ip' - ip% 
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To explain the estimates for the angle condition, we write the normal 
vector as a multilinear form on DF^: 

n+l 

N{DF) = Jn{DF) := (-1)" ^(-iy-^(£)F^A. . . £)>^A. . . DF"+^)ei e M"+^ 

i=l 

[DF^ omitted in the i*'^' term of the sum), where DF^ G for i = 1, . . . , n+ 
1 and we identify the n-multivector in with a scalar, using the standard 
volume form. The angle condition has the form: 

and we set: 

B{DF) := (3^\Jn{DF)f - {Jn{DF), en+if, 
with linearization at DFq = [In\Duo]: 

CoB[DF] = 2/3\jniDFo),DJn{DFo)[DF])-2{Jn{DFo),en+i){DJniDFo)[DF],en+i). 

Under the assumption F^ = F'^ at t = 0, we need an estimate in 
C^+^'H^ for: 

B{DF^,DF^) := B{DF^) - B{DF^) - CoB[DF^ - DF^] 

= pWJn{DF^)\^ - \Jn{DF^)\^ - 2{Jn{DFo),DJn{DFo)[DF^ - DF^])) 
-{{Jn{DF^),en+if-{Jn{DF^),en+if-2{Jn{DFo),en+i){DJniDFo)[DF^-DF%en+i)). 

It will suffice to estimate the expression in the first parenthesis; the second 
is analogous. 

We need the following algebraic observation: if Tq = [In\Duo] and T are 
n X (n + 1) matrices, the expression: 

|j„(ro + r)|2 - |j„(To)|2 - 2(j„(ro),z)j„(To)[r]) 

is a linear combination (with constant coefficients) of terms of the form: 

where the p^2) {T) are polynomials in the entries of T (with constant coeffi- 
cients), with terms of degree: 2 < deg < 2n. 
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Thus B{DF^,DF'^) is a linear combination (with constant coefficients) 
of terms: 

uoiP(^2){DF^ - DF^), uoiUojP(2){DF^ -DF^), p^2){DF' - DF''), 

with the p(2) ^ described; and hence is a hnear combination of terms of the 
form: 

(where 2 < d < 2n, 1 < j < n + 1,1 < i,l, k < n), which we write symboh- 
cahy as: 

B{DF\DF^)^ b{x){DF^ - DFY, 

2<d<2n 

where b{x) is constant or noi(x) or uoi{x)uoj (x) . For the degree d terms 
b{x){DF'^ - DF'^f, it is not hard to show the bound: 

||G^'^^||i+a < c||6||i+a||i^' - fYt"", 2<d<2n. 
We conclude: 

Lemma 7.3.Assume F, Fq, F^,F'^ are in C^+"'i+"/2(Q^; M"+i) and have 
the same initial values. Then: 

\\B{DF,DFo)\\i+^ < c(l + ||no|p)(l + ||F - Fo|p"-2)r"||F - Fo|p. 

\\B{DF\DF^)\\i+a < c(l + |ko||^)(l + \\F^ - fY''~^)T''\\F^ - 
with c depending only on Fq. In particular, if F € Bj^: 

\\B{DF,DFo)\\i+^<coR^T^, 

and if F^,F^ e Sg: 

\\BiDF\DF^)\\,+a < coT"||F^ - F\ 
with Co depending only on Fq. 
8. Local existence. 

Given a C^"*"" graph Eg over a uniformly (7^^" domain Dq C IR" (for ar- 
bitrary a G (0, 1)) satisfying the contact and angle conditions, let ipo G 
Diff^^°' be the diffeomorphism given by lemma 4.1 (with the 1-jet of 
the identity at ODq and 2-jet determined by the mean curvature of Sq at 
dDo). Then find uq G C2+"(Do) so that Fq = [(^o,wo] G C2+"(i:>o; M"+i) 
parametrizes Sq over Z)o (a = < a). 
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(Precisely, if [z,uo(^)] parametrizes Sq as a graph, and ^po is given by 
lemma 4.1, let uq = u^o 990; so uq G C^+".) 

We obtained in section 7 all the estimates needed for a fixed-point argu- 
ment in the set: 

Bl = {Fe C2+°'i+-/2(Q^,M"+i); ||F - Foil < R,F\t=^ = Fq}. 

Choose R < 1 and Tq < 1 small enough (depending only on Fq) so that, 
for F £ , F{t) = [(p{t),u{t)] defines an embedding of Dq, with ip{t) a 
diffeomorphism onto its image D{t). Let K C Imm(M", IR"+i) be a compact 
set containing DF(z) for all F G 2; G Q"^". Now consider T < Tq. 

Given F G -B^, solve the linear system (LPDE/LBC) (with initial data 
Fo) to obtain F G C2+"'i+"/2(Qr), (^h is is possible since the complemen- 
tarity and compatibility conditions hold for the linear system.) This defines 
a map ¥ : F F. 

From linear parabolic theory (e.g. [3], thm VI.21]): 

||F - Foil < M(||.F(F,Fo)|U + | |^(F»F, F»Fo)| |i+a + \\n{Dif, Dipo)\\i+a), 

where M > depends on the (7"'"/-2 norm of the coefficients of the linear 
system, that is, ultimately on ||Fo||. 

From lemmas 7.2-7.4 in section 7, it follows that: 

||F - Foil < Mco(i?T"/2 + i^^T") < R 

provided T is chosen small enough (depending only on Fo.) Thus F maps 
to itself. 

Similarly, if F(F*) = F* for i = 1,2, standard estimates for the linear 
system solved by F^ — F^ give: 

||f1-F2|| <M(||.F(F\F2)||^^||gp^i^^^2)||^^^^pp-i^^-2)||^^^) 

Again the estimates in lemmas 7.2-7.4 imply: 

_ i?2|| < Mco(r"/2 + t°)||f1 - F^ii < iiiJ?! _ i?2||^ 

assuming T is small enough (depending only on Fq). This concludes the 
argument for local existence. 

Theorem 8.1. Let Sq C M"+^ be a C^+" graph over Dq C satisfying 
the contact and angle conditions at ODq. With a = a^, there exists a 
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parametrization Fq = [ip(j,uo] G C^"'"°(Do) of Sq, a number T > depending 
only on Fq and a unique solution F G C^+"'^+"^/^((5'^; of the system: 



with initial data Fq. For each t G [0, T] , F{t) is a C^+" embedding parametriz- 
ing a surface T,t which satisfies the contact and angle conditions and moves 
by mean curvature. In addition, F(t) satisfies the orthogonality condition 
at dDo. 

The hypersurfaces are graphs. For each t G [0, T], ip{t) : Dq — D{t) 
is a diffeomorphism and = graph{w{t)), for w{t) : D{t) — > M given by 
w{t) = u{t) o if-^{t). (We have G C2+"' 'less regular' than 
or ^{t).) D{t) is a uniformly C^'^°' domain. 

Remark 8.1. This theorem does not address geometric uniqueness of the 
motion, given Sq. It only asserts uniqueness for solutions of the parametrized 
flow (including the orthogonality boundary condition) in the given regularity 
class. 

9. Rotational symmetry. In this section we record the equations for 
two rotationally symmetric instances of the problem: (i) Dq and D{t) are 
disks, and « > ('lens' case); (ii) Dq and D{t) are complements of disks in 

('exterior' case). For simplicity we restrict to n = 2. 

Let F{r) = [ip{r),u{r)\ parametrize a hypersurface S, where ip{r) = 
(j){r)er is a diffeomorphism onto its image. Here er,eg are orthonormal 
vectors, outward normal (resp. counterclockwise tangent) to the circles 
r=const. The unit upward normal vector and mean curvature are: 



{ 



dtF - g^^{DF)didjF = 0, F= [^, u] 

U\dDo = 0, N^+\Dip, Du)\QD, = P, {D^f, dnip)\dDo = 



iV = 



[-Urer,(t>r] 



1 



(0^A1(0r, Wr)[-D^'u] — {Urer,M.{^r-,Ur\D'^'f\)), 



H 



((/.2 + n2)3/2 



where: 




Simplifying: 



H = 



1 



UrCPrr + {(Pr + ^r)~rJ- 



(02+^2)3/2 
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Now consider the time-dependent case F{r,t) = [(j){r,t)er,u{r,t)]. Prom 
the above expressions, one finds easily that the equation {dfF, N) = H takes 
the form: 

~ -H^r^-^i^r,Ur)[D'^u]) = Ur{er,(pt - --T^^M{(f>r,Ur)[D'^(p]). 

In 'spht gauge', we consider the system: 

/ ^* ~ ^^^J^{'l>r,Ur)[D'^u] = 0, 

Note that 0(r, t) = r solves the ^ equation, and that in this case the u 
equation becomes: 

wt - 3 = 0. 

This can be compared with the equation for curve networks: 

^* ~ 1 , 2 = 



The boundary conditions are easily stated (we assume Dq is the unit 
disk or its complement). 

The 'contact condition' at r = 1 is tt = 0. For the 'angle condition' at 
r = 1, we find: 

Assuming 0^ > at r = 1, this resolves as: 

+ /3o<?!'r- = at r = 1 (lens case); 

(3ur — l3Q(j)r = at r = 1 (exterior case). 

(For lenses, one also has at r = 0: = and (f)r = 1.) Thus in both cases 
one can work with linear Dirichlet/Neumann-type boundary conditions. 

One reason to consider the exterior case is that (unlike the lens case) 
it admits stationary solutions. Geometrically, one just has to consider one- 
half of a catenoid, truncated at an appropriate height. For example, for 120 
degree junctions the equation for stationary solutions: 

3^ + = in {r > 1}, 
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admits the explicit solution: 

u{r) = ^(ln(2r + \/4r'^ - 3) - In 3), r > ^3/2. 



Problem. It would be interesting to consider the nonlinear dynamical 
stability of this solution (even linear stability is yet to be considered.) One 
may even work with bounded domains, by introducing a fixed boundary at 
some i? > 1, intersecting the surface orthogonally (see Section 10.) 

10. Fixed supporting hyper surfaces. Extending the local existence 
theorem to the case of hypersurfaces intersecting a fixed hypersurface S 
orthogonally presents no essential difficulty. The case of vertical support 
surface leads directly to graph evolution with a standard Neumann condition 
on a fixed boundary; we consider the complementary case where 5 is a graph. 
Let S C M"'+^ be a embedded hypersurface (not necessarily connected), 
the graph over P C M" of B G C^{V), oriented by the upward unit normal: 

Hy)-=—Hy), %) := [-i)i?(y),i] g x m, vb ■.= VT+\db{^. 

Vb 

V is assumed to be nowhere vertical in P [DB 7^ 0). To state the problem 
in the graph parametrization, we consider a time-dependent domain C 
M"" with boundary consisting of two components diD{t) and d2D{t), both 
moving. The hypersurface is the graph of w(-,t) over D[t), solving the 
parabolic equation: 

wt - g'\Dw)wij = in £ := (J D{t) x {t] e M"+^ x [0, T], 

te[o,T] 

with boundary conditions: 

w{-,t)\Q^D{t) = 0, ^/l + \Dw\'^\Q^JJ^^^ = 1/(5 
(as before), and on (?2-D(i): 

w = B, Vw VB = -1. 
(The first-order condition on d2D(t) is equivalent to {iy,N) = 0). 

Differentiating in t the boundary condition w = B leads easily to an 
equation for the normal velocity of the interface T{t) = d2D{t): 

. vH 

J- n — J-, 

Bn - Wn 
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Note that Wn at d2D{t) can be computed from Bn-, since: 

-1 = • = WnBn + |V^5p; 

in particular neither Bn nor Wn can vanish (so both have constant sign 
on connected components of ^2-0), and one easily computes: Wn — Bn = 

-V%/Bn. 

Let A = S n 5 be the (n — l)-manifold of intersection, the graph of w 
(or B) over d2D. Given the graph parametrizations of S and S: 

G{y) = [y,w{y)], M{y) = [y, B{y)], y G 

and T G Td2D, we have the tangent vectors: 

Gn ■■= [n, Wn] G TS, Gb := [VS, -1] = -vbi^ G TS, := [r, V-uj-r] G TA, 
and the second fundamental forms of S and 5 (for e G M" arbitrary): 

A(dGe, dGe) = -d'^w(e, e), ^(dBe, cflBe) = —(fB{e, e). 

V VB 

From {i^jN) = at d2D, it follows easily that (cp. [13]): 

A{Gr, v) = -A{Gr, N), T G TdD. 

For the remainder of this section, we concentrate on the boundary con- 
ditions at To establish short-time existence, we consider as before the 
parametrized flow: 

Ft-trgd^F = 0, g = g{dF) , F = [99, n] . 

The contact and angle boundary conditions are: 

Again we have two scalar boundary conditions for n + 1 components. Here 
the solution is easier than at the junction. With the notation Fn = dFn = 
[(PmUn\, we replace the angle condition by the 'vector Neumann condition': 

Fn -L T5, or F„ = -avsi^ on d2DQ, 

where a : 92-Do — > K, or equivalently (since this leads to q = —Un)'- 

= -Un{VB O If) on 92-Do- 
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Clearly the Neumann condition implies the angle condition {N, v o ip) = 

0, but not conversely. This linear Neumann-type condition can easily be 
incorporated into the fixed-point existence scheme described earlier. 

There is one issue to consider: the and Ist-ordcr compatibility con- 
ditions must hold at at t = 0. The initial hypersurface Sq uniquely 
determines wq and Dq C (satisfying wq = B and Vwq ■ VB = — 1 on 
52-Do)) and then once (^o € Diff^Do) is fixed, uq = wo o (pQ is also deter- 
mined. We may assume: 

(fo = id, (fion = V-B on 52-Do, 

so: 

UOn = ■ VOn = Vwq • V-B = -1 on d2Do, 

and then the Neumann condition Fgn^g^Oo ~ holds at t = 0, on 820^. 

The first-order compatibility condition is: 

trgcfuQ = ut = V-B • (fit = V-B • trg(f(pQ on ODq, 

or equivalently: 

trg{i^,d^Fo) = on ODq. 

(This is not a mean curvature condition; the mean curvature of T,q is H = 
trg{N,d''Fo).) 

From now on we omit the subscript 0, but continue to discuss compatibil- 
ity at t = 0. First observe that the Neumann condition leads to a splitting of 
the induced metric. Given r G Td2DQ, let Ft = dFr G TA. Then (recalling 
Un = -1 on 52-Do): 

{Fr, Fn) = ([r, dBr], [y^n, Un]) =VB-T-VB-T = 0. 

Thus we have: 

trg{,y, d^F) = g''\v, d^F{Ta, n)) + ^""(j., d^F{Fn, Fn)), 
for a local basis {Ta = dFTaYaZl of TA, with gab = {Ta,Th) and gnn = 

\F |2 - 7;2 

Differentiating in n the condition Un = ■ (pn ( assuming, as usual, 
n extended to a tubular neighborhood M of 52-Do as a self-parallel vector 
field), we find: 

Unn = cPw{n, V-B) -|- Vw ■ d^(p{n, n). 
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This is used to compute: 

{u,(i^F{n,n)) = — [■u„„ - VB ■ (f(p{n,n)] 

VB 

= —[(fw{n, VB) + {Vw - VB) ■ (f(p{n, n)] 

VB 

= -vA{Gn, 1^) + —{wn - Bn)n ■ (f(f{n, n). 

VB 

Bearing in mind the expression for Wn — B^ found eaxher, the compatibihty 
condition may be stated in the form: 

• d\{n,n) = -vA{Gn,u) + g''\d^F{Ta,n),i^). 

We are now in the same situation as in section 4: given the 1-jet of <fo 
on d2Do, we extend (po to a tubular neighborhood M of 92-Do (and then 
to all of Do), so that n ■ d!^ip{n,n) has on 92-Do the value dictated by the 
compatibility condition (using Lemma 4.1(ii)). We just need to verify that 
the right-hand side of the above expression depends only on Sq, S, and the 
1-jet of (fio over 52-Do- Clearly only the term ^"^(i/, (PF{Ta, r^)) is potentially 
an issue. 

Fix p G Sind let {Ta} be an orthonormal frame for Tc?2-Do near p, 

parallel at p for the connection induced on 92-Do from M". If /C denotes the 
second fundamental form of 92-Do in I^"; we have: 

Tain) = JC{Ta, n)n ( at p) 

(on the left-hand-side, Tt, is regarded as a vector- valued function in M"). Still 
computing at p, this implies: 

(fF{ra,n) = TaidFn) - dF{Ta{n)) 

= Ta{dMn) - /C(ra, Tb)Fn 
= d^M{Ta, n) + lC{Ta, T;,)]B„ - /C(Ta, T„)F„, 

where F„ = —vu and B„ = dMn G TS. Hence: 

{u, d^Fija, n)) = {ly, d^MiTa, n)) + v/C(r„, n) = A{Ta, n) + v/C(r„, n). 
This clearly depends only on 5 and on Sq. We summarize the discussion in 
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a lemma. 

Lemma 10.1 Let Sq = graph{wo) be a graph over Dq C (a 
uniformly domain), intersecting a fixed hypersurface S = graph{B) over 
ODq. Consider the parametrized mean curvature motion with Neumann 
boundary condition: 

F eC^'^{Dox[0,T]) F=[ip,u] 

Ft - tVgSF = 0, g = g{dF), uo ip = B and Fn -L TS on dD^. 

Then ipQ G Diff{Do) can be chosen so that (with uq = wq o ipo) the initial 
data Fq = [ipo , uq] satisfies the order zero and the first-order compatibility 
conditions at t = and ODq: 

ipOn = -UoniVB O ipq), {v O ifQ, trg^cfFo) = 0. 

Remark 10.1. Differentiating dwTa = dBTa along T(,, we find: 

d^w{Ta, n) - d^B{Ta, u) = {Wn - -B„)/C(ra, u) 

(reminding us that, although it; = B on dDo, the tangential components of 
their Hessians do not coincide.) Prom this follows the expression for /C in 
terms of A and A: 

K^{Ta,n) = ^—B-b'MTa,Tb) - VBA(Ta,Tb)]. 

Wn - Bn 

It is also easy to express the corresponding traces in terms of the mean 
curvatures and Ti^ of A in E and S: 
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11. Boundary conditions for the second fundamental form. 

To understand the long-term behavior of a graph (S^) in M""*"^ moving 
by mean curvature and intersecting R" at a constant angle, we need to 
consider the evolution of its second fundamental form. Working in the graph 
parametrization the boundary conditions are easy to state (and linear): 

n A ^0 

W\dD{t) = U, anW\QD(t) = 

where n = nj is the inner unit normal to dD{t). It is possible to reparametrize 
the Tit over a different time-dependent domain 'D{t), obtaining mean curva- 
ture flow: 

rt-'^{t)^W+^, dtJ' = HN, 
with boundary conditions: 

For this parametrization the evolution equation for the second fundamental 
form (and its covariant derivatives of arbitrary order) is well-understood [5]. 
The disadvantage is that the unit normal N\jr^ depends non-linearly on the 
components of and as a result the boundary conditions for the second 
fundamental form (which are needed for global estimates over spacetime 
domains) do not admit simple expressions. Therefore we choose to work 
with graph flow, at the cost of having to derive and understand a new set 
of evolution equations. The equations for h and the mean curvature H are 
derived in an appendix. In this section we derive boundary conditions; the 
development is similar to work of A. Stahl [13] for MCF of hypersurfaces 
intersecting a fixed boundary orthogonally. 

It is easy to see that h splits on dD{t): if r G TdD{t) is a tangential 
vector field, and n = nt is the inner unit normal: 

hin, r) = -cfw(n, r) = -(T(wn) — Dw ■ Vrra) = on dDit), 

V V 

since = /?o//3 on the boundary and V^n G TdD{t) (V is the euclidean 
connection.) In particular, it follows that h{Dw,T) = on dD{t). 

Remark. Already this simple fact cannot be shown for a(z^, r), the second 
fundamental form in the MCF parametrization, regarded as a quadratic form 
on P(i). 
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Boundary condition for H. We derived in section 2 the equation for the 
normal velocity of the moving boundary Tf = dD(t): 

tn = -—H = -^H aidDit). 

Since (AT, e„+i)(r(t), t) = ^ on dD{t) we have: 

{dtN,en+i) = -{dkN,en+i)t^ 
where dkN = —g^^hikGj, with e„+i component: 

{dkN, e„+i) = -g'^Wjhik = -^h{Dw,dk) = -^Wnh{n,dk). 
Hence we find, on dD{t): 

{dtN, en+i) = ^h{n, f ) = ^tnh{n, n) = -pHhnn- (11.1) 

(We set hnn '■= h{n,n)). Denote by the gradient of Sj, in the induced 
metric (V^/ = g^^fiGj). Using dtN = —V^H — Hv~^V^v, combined with 
the expressions (valid on dD{t)): 

, E , _ VnWn _ ^1 _ 2, 

yZ 

we find on dD{t): 

{dtN, en+l) = -P(3o{Hn + PoHhnn). (11-2) 

Comparing expressions for (11.1) and (11.2) for {dtN, e„+i) yields a Neumann- 
type condition for H. We state this as a lemma (including the evolution 
equation derived in the appendix). Here L = Lg denotes the operator 
L[f] = dtf — trgD^f and lo = Dw/v, a vector field in D{t). 

Lemma 11.1. For the surfaces evolving by graph mean curvature 
motion with constant contact angle, the mean curvature satisfies: 

L[H] = \h\lH + Hh'^{LU, Lj) - H^h{u;, lo) on D{t) 

dnH = ^Hhnn On dD{t) 



28 



Boundary conditions for hij. Fix p G dD{t) and let (tq) be an orthonor- 
mal frame for TpdD{t) (in the induced metric), satisfying V^^Tf,{p) = (V^ 
is the connection induced on Tf by the euclidean connection d, or, equiva- 
lently, by V, the Levi-Civita connection of the metric g in D(t)); we extend 
the Ta to a tubular neighborhood of Tf so that V^Ta = 0. Differentiating 
h{n, Tb) = along Ta, we find: 

{VTah){n,Tb) = -/i(V^„n,rb) - /i(n, V^^Tf,). (11.3) 

The second fundamental form /C(r, r') of Ft in (D(t),eucl) (equivalently, in 
{D(t),g)) is defined by: 

draU = V^^Tfe + JC{Ta, n)n On dD{t). 

To relate /C to h\dD(t): "^ote that since u; = on dD{t): 

hiTa,n) = ^d'^'w{Ta,n) = ^{Ta{T},w)-Dwdr^T},) = -dr^Tf,-^ = - PolC{Ta, u) . 

(So we see that convex with respect to n corresponds to concave over 
D{t), as expected). In the appendix we observe that Vg.Sj = {hij/v)Dw. 
Then: 

Vr„r6 = Ti{{Tl)idj + TlVQ^dj) = dr^n + -TlTlhijDw 

= ^\a'^h+ K^{Ta,Tb)n+—h{Ta,Ti,)n = {- ^ + Po)h{Ta, n)n = -—h{Ta,n)n 

V po po 

at p, given our assumption V^^Ti,{p) = 0. We use this immediately to 
compute, at p: 



13^ 1 
^Ta^ = {^Tan,n)gTb = -{n,Vran)gn = -;r\n\lh{Ta,Tb)Tb = --rh{Ta,n)Tb, 

Po Po 



since |n|^ = gijU^n^ = 1 + w^ = (3 ^ at p. Using these expressions for Vr„n 
and Vr^Tb in (11.3) above and recalling the Codazzi equations: 

1 p'^ 

{Vnh){Ta,Tb) = {V r^iu^Tb) = - -5- V /i(ra, Tc)/i(Tc, Tfo) + — h{Ta , n) hnn ■ 

PO ^ Po 
This can also be written in the form: 

Po(ynh){T,T') = -{h'n\r,T')+P^KnhiT,T'). (11.4) 
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It turns out the expression for the n-directional derivative of /i(r, r') is ex- 
actly the same (at dD{t)): 

MniKry)) = -{h"^''f{T,T')+(3'hnnh{T,T'). (11.5) 

The reason is that V„ra = at the boundary, also for the (7-connection: 

VnTa = dn{Ta) + nWiVo^dj = + ^/l(n, Ta)Dw = 0, 

so in fact: 

(V„/l)(Ta,rb) = n{h{Ta,n)) = dn{h{Ta,n)). 

As done in [13j, we combine this with the result for Hn to compute 
{SJ nh){n,n). From: 

Hn = Vnitrgh) = trgiVnh) = p\Vnh){n,n) + ^(V„/i)(T„, r„). 

a 

Here we used \n\g = on dD{t), which also implies H = fPhnn + 
T.aKra,Ta). Using also = E(/i*"")'(^a,ra), we find for (V„/i)(n,n): 
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(3^ 
Po 


= ^(|/i*'^"|2 

/3o 


+ /3'0 = 


Po 



since y"" = /3 at dD{t). Equivalently: 



/3o(V„/i)(n,n) = on aZ)(t). 



It is easy to obtain the corresponding expression for the euclidean con- 
nection. Noting that at dD{t): 

VnU = dnU + n^n^-hijDw = Pohnnn, 

V 

we find: 

{dnh){n,n) = n{hnn) = iVnh){n,n) -F2/i(V„n,n) = {Vnh)in,n) + 2(3ohl^, 
so that: 

Podn{Kn,n)) = ^\h\l + 2/?2C on ODit). 
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We record these results as a lemma, including also the evolution equa- 
tions derived in the appendix. 

Lemma 11.2. Under graph mean curvature motion with constant con- 
tact angle, the second fundamental form satisfies the following tensorial 
evolution equations. {Cij and Cij are symmetric 2-tensors cubic in h, given 
explicitly in the appendix.) 

(i) For the operator L = Lgi 

L[hij] = -2[hUw{hjk) + h!}d^{hik)\ + Cij on D{t), 
with boundary conditions on dD{t): 

' h{n,T) = 

< (3odn{h{T,T')) = -(/i*"")2(r,r') +/32/inn/i(r,r') 
^ Podn{h{n,n)) = ^\h\l + 2p^hl„ 

(ii) For the operator df — Ag, where Ag is the Laplace-Beltrami operator 
of g: 

{dt - Ag)[h]ij = H{\7^h)ij + Hih{uj, dj) + Hjh{uj, di) + Cij on D{t), 
with boundary conditions on dD(t): 

h{n,T) = 

/3o(Vn/i)(r,r') = -{h'-ny{r,T') + P^hnnh{T,r') 
PoiV nh){n, n) = ^\h\l 

It is also useful to compute the boundary condition for \h\^. Using 
Lemma 11.2(ii), we have at dD{t): 

iPo/2)dn\h\l = Po{Vnh,h)g 

= (3oP'^{Vnh){n,n)hnn + PoY^b c(ynh)iTa,n)h{Ta,n) 

= f5^\h\lKn + Eab[-ih"'nHra,n) + P^Knh{Ta,n)]hiTa,n) 

= PWh\l + \h"'^\l)hnn-trg{h"'^f 

Since on dD{t): trgh^ = (3^ {hnnY +tr g{h^°''^Y , we may state this in a slightly 
different form. Including also the evolution equation for \h\'^ (see appendix), 
we have the following lemma. 
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Lemma 11.3. Under graph mean curvature flow, the function \h\g 
satisfies the evolution equation and Neumann boundary condition: 

r {dt - Ag)\h\l = -2\Vh\l + Hd^\h\l + 2\h\l - AHh\uj,uj) - 2H\h\lh{u^ , u:) 
\ {l3o/2)dn\h\l = 2p^\h\lKn - trg{h^) on dD{t) 



12. A maximum principle for symmetric 2-tensors. 

By the local existence theorem, for suitable initial data we have a mean 
curvature motion F = [ip, u] G C^+^'^+^/^^g^^ ]^n+i)^ vjhexe Qq = DqX [0, T] 
and, for each t G [0,T], ipt : Dq D{t) is a C^"*"" diffeomorphism. In 
particular, with 6 = a'^ , wt = ut o ip^^ : L'(t) — > M defines a graph m.c.m. 
w G C'^^^'^^^^'^ {E;M) in an open spacetime domain 

E= [j D{t) X {t} C M" X M. 
te(o,T) 

We have a (72+a,i+Q!/2 diffeomorphism: 

$:Qo^^, Hx,t) = {ipt{x),t), 
which, for any > 0, restricts to a diffeomorphism Qt^ Et^, where: 
Qt, =DoX (to,T), Et, = U D{t) X {t}. 

The parabolic boundary of E is the disjoint union of 'base' and 'lateral 
boundary': 

dpE = {Do X {0}) U diE, diE = U dD{t) x {t}. 

te{o,T) 

(The notions of 'parabolic boundary, 'base' and 'lateral boundary' have gen- 
eral definitions for arbitrary bounded spacetime domains (see ^), but using 
$ it is easy to see that they are given by the above sets.) In particular, note 
that $ defines a diffeomorphism (C''^"'""'^^"'"")/^) up to the lateral boundary, 
if Do is a C7'=+° domain and F G (:7*=+"'(*^+°)/2)(Qo)): 

Qto U diQto Etg U diEtg, 

for each to > 0. 
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Denote by L the operator L = dt — g^^ {Dw)didj , so Lw = in E and 
w = on diE. The fohowing height bound is immediate. 

Lemma 12.1. Assume < wq < M in Dq. Then {) <w < M in E (and 
vanishes only on diE). 

Proof. Follows from the weak maximum principle for the operator L, 
since < u; < M holds on the parabolic boundary dpE. 

It is well-known that the function v = ^/l + \Dw\^ solves the evolution 
equation (assuming Dw G C^'^(£'), see e.g. ^): 

2 2 

L[v\ + -g'^ViVj = -v\h\l, or L[v\ = — \Dv\l - v\h\l. 

From the maximum principle, we have the following global bound on v 
(equivalently, on |-Dii^|). 

Lemma 12.2 Assume w is a solution with Dw G C'^'^{E). Then we 
have on ^: ^ 

v{z) < max{supD(^to)v{x,to), -}. 

Proof. By the weak maximum principle, max^ v = maxg^E v. Note 

It follows from this lemma that gij{t) is uniformly equivalent to the 
euclidean metric in D(t): if v < v in E, and X is a vector field in D{t): 

\X\l < \X\l = g^jX'X^ = \X\l + {X • Dwf < \X\1{1 + \Dw\^) < v^\X\l. 

Also, if uj := v~^Dw: 

. ,2 \Dw\l ^ 1 ^ 1 

yZ yZ yZ 

The main result in this section is a maximum principle for symmetric 
2-tensors satisfying a parabolic equation on a spacetime domain such as E 
(image of a cylinder under a diffeomorphism of the special type $). 

Recall the 'boundary point lemma' for scalar equations . It holds for 
open spacetime domains Q C x satisfying an 'interior ball condition': 

Interior ball condition: For each P = (p, t) G difl we may find a ball 
B (in the euclidean metric in M""''^) which is tangent to diQ, only at P and 
satisfies: (i) the line segment from P to the center of the ball is not parallel 
to the t axis; (ii) B n {t < t} C ^In {t < t}. 
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Remark. For the domain of interest the interior ball condition follows 
from the fact that diE = <^{dDo x (0,r)), with ^> G C^'^iDo x (0,r)) of the 
special form above. 

Lemma 12. 3. ([H] Thm.6, p. 174.) Let Q C M" x M+ be a connected 
open set satisfying the interior ball condition. Assume / S C^'^(ri) satisfies 
the uniformly parabolic inequality: 

dtf - trg(ff - dxf < 0. 

Here g = gtisa Riemannian metric in each section r2(t), and Xt is a bounded 
vector field in Q(t). Denote by n = nj the inner unit normal of i}(t). 

Assume the supremum M of / in ilj := 17 n {t < t} is attained at the 
point P e dQ{t), and that / < M for t < t. Then (i„/(P) < 0. 

We now state the hypotheses of our tensorial maximum principle. 

£; C X [0, T] is the image of a cylinder Dq x (0, T) under a C^'^ 
diffeomorphism <I> of the form ^{x, t) = {ipt{x),t), with tpt '■ Do D{t) a 
diffeomorphism up to the boundary, for each t E [0, T] {D{t) is the t = const. 
section of E); Dq C M" is assumed to be the image of the closed unit ball 
under a diffeomorphism.) In particular, the lateral boundary diE is of 
class C^'^. On diE we have the inner unit normal n = nt £ M". Extend 
nt to a vector field in all of D{t) so that it is in C^'^(£', M"), arbitrarily 
except for the requirements that \n\ < 1 pointwise and d„n = in a tubular 
neighborhood of dD{t) (equivalently, n^din^ = for each j.) Fix i? > so 
that D{t) C Sj?(0), for each t G [0,r] 

The assumptions on the coefficients are given next. 

g = gt '\s a, t-dependent Riemannian metric in D[t), uniformly equivalent 
to the euclidean metric for t G [0, T]; 

X = is a bounded t-dependent vector field in D{t); 

q = q{z, m) assigns to each z £ E and each m in S (the space of quadratic 
forms in M") a quadratic form q £ §. q is assumed to be C^'^ in z, locally 
Lipschitz in m (uniformly in z £ E); 

b = b{z, m) G S is defined for z £ diE, with the same regularity assump- 
tions as q. 

Remark. We state the following theorem in terms of the Laplace-Beltrami 
heat operator dt — Ag and the g- Riemannian connection V, but the result 
also holds for L and the 'euclidean connection' d. 

Theorem 12.1. Assume m £ C^'^(£^;S) satisfies in E the tensorial 
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differential inequality: 

dtruij - {^gm)ij < {Vxm)ij + qij{-,m{-)), 
and on diE the boundary condition: 

(Vnm)ij(z) > bij{z,m{z)). 

Suppose the fimetions q and b satisfy the following 'null eigenvector condi- 
tions': if, for some m G §, V G is a null eigenvector of m (TTiijV^ = OVi), 
then, for any z e E (resp. any z G diE): 

qij{z,7n)V'V^ < (resp. bij{z,m)V'V^ > 0). 

Then weak concavity of m at t = is preserved: 

m < in D{0) ^ m < in ^. 

Proof. The assumptions imply there is X > (depending only on E and 
on the functions X, g, n, q and b) satisfying: 

\n\c^AiE)<K, \X{z)\eucl < K, \g{z)\ + \g-\z)\<K, z e E; 
and if m, m G C'^'^{E, S) satisfy (for some n : E ^ M+): 

—li{z)g < m{z) — m{z) < n{z)g 
(where the inequality of quadratic forms has the usual meaning), then also: 
q{z,m{z)) < q{z,m{z)) + KiJ,{z)g, z e E, 
b{z,m{z)) > b{z,m{z)) — KiJ,{z)g, z G diE. 

Now define, for 2; G ^, z = {x,t): 

ip{z) := -2Kn{z) ■ x := 2Ks{z), 

where we use the euclidean inner product and, on diE, s is the 'support 
function' of dD{t) (positive if D{t) is convex and contains the origin). It is 
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clear we may find M = M{R, K) > depending only on K, R and |n|(72,i 
so that: 

\ip\c2,i < M, \dip\l + \Agip\ <M, \X ■ difl < M. 

We assume also M > K. Now, given m as in the statement of the theorem 
and given constants e > 0, 7 > and 5 > 0, define for := E {t < 5}: 

m{z) := m{z) - {et + ^e'^^'^)g, z e eK 

Clearly m G C^^^ [E^ ■,^) . We now derive the constraints on (5, e and 7. It 
will turn out that 8 must be taken small enough (depending only on R), 
e > is arbitrary and 7 is e times a constant depending only on K, R. 

The following inequalities are easily derived: 

q{z, m{z)) < q{z, ih{z)) + K{et + 76'^^'^))^?; 

Vxm = Vxrh + ^{e.'^dx^)g < Vxrh + {ie'^M)g; 
dtjfi = dtm — eg — {'ye^dtip)g < dtm + {'ye^M)g — eg; 
Agrh = Agd^m - -fe^iW]] + Ag(p)g > AgUi - {^e'^M)g. 
b{z, m{z)) > b{z, m{z)) — K{et + ')e^)g 

We use this to compute: 

dtrh — Agrh < dtm — Ag-m + {2je'^M)g — eg 

< q{z, m{z)) + Vxm + {2je'^M)g - eg 
< q{z, m{z)) + Vxm + K{et + ^e'^)g + (3M7e^)5 - eg 
< q{z, m{z)) + Vxm + Metg + AM^e'^g - eg 
(using < M in the last step). We conclude the inequality: 

dtifi — AgiJi < q{z, m(z)) + Vxm — {e/2)g (13.1) 

will hold in E^, provided the constants are selected so that, for z E E^: 

4M7e'^(^) + Met < e/2. (A) 

Turning to boundary points z = {x, t) G diE, note that dn(p = —2K, so 
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that: 

V„m(z) = V„m(z) - {je^^'Unip{z))g > b{z,m{z)) - {-fe'^^'Unip{z))g 

> b{z,m{z)) - K{et + je^^'^)g - {-ie^^'Un^{z))g 
> b{z,m{z)) +K{je'^^^'> - et)g, 
implying the inequahty: 

V„m(z) > b{z,m), z G diE^ (13.2) 

will hold provided the constants are chosen so that, on diE^: 

ei<7e^(^). {B) 

Bearing in mind that, on E: e~^^^ < e'^^^^ < e^^^, it is not hard to arrange 
for (A) and (B) to hold, or equivalently, for: 

Given e > 0, define 7 so that 10M7e^^^ = e. Then the second inequality 
holds, and so will the first, provided: 

et < 7e-2^^ = (e/10M)e-^^^, 

which is true for any e > 0, if 5 is defined via 6 := e~^^^ /lOM (recall 
te[0,6]). 

Note that, since m < at f = 0, it follows that m is negative-definite at 
t = 0, and hence also for small time, and we claim that this persists through- 
out E^ , so that (letting e ^ 0) m < in £'^. Restarting the argument at 
t = 6, wc see this is enough to prove the theorem. 

To prove this claim, suppose (by contradiction) m acquires a null eigen- 
vector / y G M" at a point zi = {xi,ti) G E^ , with ti £ (0,(5] the first 
time this happens. 

Let f{z) := rhijV^V^, z e E^ (that is, we 'extend' V to E^ as & 
constant vector.) It follows from (13.1) that / satisfies in E^: 

dtf < {Agm)ijV'V^ + {Vxm)ijV'V^ + qij{;m)V'V^ - 

A short and standard Riemannian calculation (using the fact that F is a 
null eigenvector for m) shows that: 

dxf = {Vxm)ijV'V^, Agf = {Agm)ijV'V^. 
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Using the null eigenvector condition for q, we find that / satisfies in the 
strict inequality: 

dj <trgcpf + dxf. 

This shows xi cannot be an interior point of D{ti), for then (as a first-time 
interior maximum point for /) we would have Agf{zi) < and df{zi) = 
0, contradicting dtf{zi) > 0. Thus xi G dD{ti). Since / satisfies the 
differential inequality just stated and zi = {xi^ti) is a first-time boundary 
maximum m £"5, the parabolic Hopf lemma (Lemma 12.3 above) implies 
dnf{zi) < 0. On the other hand, as seen above in (13.2): 

dnf ={ynrh)ijV'Vi > hij{zi,m{^i))V'V^ > 0, 

from the boundary null-eigenvector condition. This contradiction concludes 
the proof. 

Corollary 12.2. Suppose m G C'^'^{E,S) satisfies the same differential 

inequality, with the same hypotheses on the coefficients as in the theorem 
(including the null eigenvector condition for q), and the boundary conditions: 

m{z){n, r) =0, Vz = (x, t) G diE, r G T^dD{t) 

{Vnm){n,n) > bnn{z,m{z)) 

(V„m)(T, r) > h'''^{z, m{z)){T, r), r G T^dD{t), 

for functions bnniz,rh) from diE x S to M and 6*"" assigning to (z, m) G 
diE xE>, z = {x,t), a quadratic form in TxdD{t). Suppose hnn > in £J x S 
and 6*"" satisfies for each rh G S: 

rhijT' = for some r G T^dD{t) =^ b^"'''{z,m){T,T) > 0. 

Then, as in the theorem, weak concavity is preserved: 

m<0ati = 0^m<0in^. 

Proof. This is proved as the theorem, with the following change in the 
last part of the proof: if / F G M" is a null eigenvector of m (defined as 
in the proof of the theorem) at a boundary point zi = {xi,ti) G diE, write: 

V = V''n + V'^, eT^,dD{ti). 

Assume first ^ 0. Then (noting that rh splits at the boundary if 
m does), we see that ra is a null eigenvector of m at zi, so we define 
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f{z) = mij{z)n^ {zi)n^ (zi) and repeat the argument. At zi, (Vn?Ti)(n, n) > 
bnnizi,rh{zi)) > leads to a contradiction with the parabolic Hopf lemma, 
as before. 

If = 0, then V'^ S Tx^dD{ti) must be a null eigenvector of m 
at the boundary point zi, and then we run the argument with f{z) = 
m{z){y'^ ,V'^)^ leading to a contradiction, as before. 

Corollary 12.3. Let w G C'^''^{E) define a MCM of graphs with 
constant-angle boundary conditions, where E is as in the statement of the- 
orem 12.1. Then weak concavity is preserved: 

/i<0att = 0^/i<0in£;. 

Proof. From Lemma 11.2, h satisfies: 

{dt - Ag)hij = HVujhij + q{z, h)ij, 

q{z, h)ij = Hih{uj, dj)+Hjh{iv,di)+\h\'^ghij+Hh{di,uj)h{dj,uj)-Hh{ui,ui)hij. 

(Here Hi, Hj, H and are regarded as fixed functions oi z £ E.) Clearly q 
satisfies the null eigenvector condition, since qijV^V^ = when hijV^ = 
for all i. In addition, expressions obtained for in lemma 11.2 show that 
the boundary conditions in Corollary 12.2 are satisfied, with: 

bnniz,m) ^ 0, = -((rn)*'^")^ + /J^^^^m*'^-. 

Hence the claim follows from Corollary 12.2. 

Remark 12.1. For less regular solutions, we may apply the theorem to 
a domain Et^^ = E r\ {t > to\, for arbitrarily small 5 > 0. Thus, assuming 
/i < at t = (strictly negative-definite), we conclude from corollary 13.3 
that /i < for all t. 

Remark 12.2. It seems plausible that a slightly different version of the 
result in this section could be used to strengthen the conclusions of [13]. 

Finite existence time. 

It is not difficult to derive that the flow is defined only for finite time in 
the concave case. 

Lemma 12.4. Let w E C'^'^{E),E C M" x [0,r) define a graph MCM 
Tjt with constant-angle boundary conditions on a moving boundary. Assume 
So (and hence Sj, for all t) is weakly concave. Then: 
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Assume -ff|t=o < -f^^o < (where Hq is a negative constant). Then 
T <U = ^TTTT— ( we are assuming T = sup{t e [0,r);D(i) ^ 0}). Here 
c„ > depends only on n and an upper bound for v in 

The proof is based on the evolution equation and boundary condition 
for H (see Appendix 2: lo = Dw/v): 

L[H] = \h\lH + Hh^{u,u) - H^h{u;,u), = {P^/po)Hhn„. 

Since h?{oj,oj) > 0, \h\l > {l/n)H'^ and (given that h < 0) h{uj,uj) > 
\Dw\^H, we have: 

L[H] < ^H^ + \Dw\'^H^ < CnH^, 

where Cn depends on n and on sup^; \v\ (already known to be finite). Let 
0(t) solve the o.d.e. = Cn(l)^,<^{0) = Hq: 

cj){t) = Ho[l - 2cnHlt]-^'^, 0<t<U:= — 

Then with := {l/n){H'^ + H(p + 4?) > 0, setting x = H-cf): 

L[x] < V'X in E; 

Xn = -rix + 4>)hnn > -^X On diE 
PO PO 

(since (f) < and /7,„„ < 0). Given that x < at t = 0, it follows from the 
maximum principle that % < 0, or < ^ in [0, min{T, t^}). This shows 
< T is impossible, since cf) ^ —oo as t ^ t^. 

Remark 12.1. It would be natural to try to show that a negative upper 
bound Hq on the mean curvature (at t = 0) is preserved, at least under 
the assumption of concavity. Unfortunately, the evolution equation for H 
(under graph m.c.m.) does not lend itself to a maximum principle argument. 
Letting u := H — Hq, we have 

L[u] = \h\lu + uh'^{uj,uj) - u{H + Ho)h{LO,Lo) + HqQ in E, 

Q := \h\l + h'^{u,io) - HQh{uj,uj). 

At a point where n = 0, we would need to show L[u\ < 0. But it is not true 
that Q > at such a point. («„ > does hold at boundary points.) 
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The exception is if n = 2 (under an additional condition). Let Co = 
oj/\oj\g,Cj = oj-^ /\oj\g. It is easy to check that B = {a),tD} is a g-orthonormal 
frame at each point where u) ^ {). Then with: 

a := h^LO , Lo) , b := h^LO , Co) , c := h{uj , uj) , 

we have: 

/i^(a),w) — Hh{u,Lb) = + 6^ — (a + c)a = b'^ — ac = —A, 

where A, the determinant of the matrix of h in B, is non-negative if /i < 0. 
In particular: 

h^{uj,uj) - Hh{uj,uj) = -l^^lgA < 

in the concave case. Now consider the expression denoted above by Q (at a 
point where n = 0, or = Hq). Since \uj\^ = \Dw\^: 

Q = \h\l + h^{Lo,Lo) - Hh{uj,uj) = a^ + 2b^ + + \Dw\^{b^ - ac) 

= 6^(2 + \Dw\'^) + - \Dw\'^ac + c^, 

so Q > provided iDwp < 2. This last condition is equivalent to v < \/3, 
and hence (Lemma 12.2) is preserved by the evolution if it holds at t = 0. 
Thus we have: 

Proposition 12.5. Assume n = 2, h < and v < \/3 on Sq (in 

particular, /3 > l/Vs). Then H < Hq < at t = implies H < Hq for all 

t S [0 , Tqjiax ) • 

13. Global bounds from boundary bounds for V"/i. 

In this section we begin to develop a continuation criterion for solutions 
of graph mean curvature motion with constant contact angle, based on the 
second fundamental form. Our first observation is that the supremum of \ h\g 
on the moving boundary controls its value in the interior. Recall we already 
have a bound on sup^;^ (Lemma 12.2); it is a well known-fact for mean 
curvature flow of graphs that this implies interior bounds for the second 
fundamental form and its covariant derivatives ([Ij, ^). In the next lemma 
we describe a global bound for mean curvature motion of graphs with moving 
boundaries. 

Lemma 13.1. Let w : — > M be a (sufficiently regular) solution of 
graph m.c.m in a spacetime domain E C M" x [0, T], where T < oo. Assume 
the first derivative bound v{x, t) < v holds globally in E. Then if the bound 
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\h\g < ho holds on the paraboUc boundary dpE, we also have the global 
bound: 

\h\g < ao in E, 

for a constant oq depending only on n,v, hQ,T and the initial data of w. 

Proof in the concave case. The proof is simpler under the assumption 
that h is negative-definite. (As shown in the previous section, this is pre- 
served if it holds at t = 0). 

Notation: In this proof, the norms of tensors in D{t) are always taken 
with respect to the induced metric g, so we write e.g. \h\ for \h\g, |V/i| for 
\Vh\g, and = g'^ fifj for a function /. 

Recall the evolution equations: 

L[v] = -v\h\^ - -\Dv\'^, so L[v^] = -2v^\h\'^ - 6|L>uP, 

V 

L[|/i|2] = -2|V/ip + 2|/i|^ - 4Hh^{u;,u;) - 2H\h\^h{uj,uj). 
In the concave case, since H < and h^ is negative-definite, this implies: 

L[\h\^] < -2|V/i|2 + 2|/i|^ 

The idea of the proof is to apply the maximum principle to / = In 
the evolution equation for /: 

L[f] = v'L[\hf] + \hfL[v'] - 2{D\hf,Dv')g, 

the terms ib2t;^|/i|^ cancel exactly, and we have the inequality: 

L[f] < -2v'^\Vh\'^ - 6|/i|2|Dt;|2 - 2{D\h\'^ , Dv'^)g. 

The term with the inner product can be estimated in two ways: 

\{D\h\^,Dv'^)g\ < \D\h\^\\Dv^\ < 4\h\v\Vh\\Dv\ < 2v^\Vh\^ + 2\h\^\Dv\'^ 

and: 

{D\h\\Dv^)g = -{Di\hfv'),Dv')g~^-^\Dv^\^ = -^{Df,Dv^)g-m'\Dv 
Using the second expression, we have: 

L[f] < -2v^\Vhf-6\h\^\Dv\^-^{Df,Dv^)g + 4\h\^\Dv\^-{D\hf,Dv^)g, 
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and then estimating the remaining inner product term from the first expres- 
sion: 

L[f] < -2t;2|V/i|2-6|/ip|Z)wp-i^(L>/,L'?;2\ ^4|/^|2|^^|2^2t;2|V/i|2+2|/i|2|I)t;|2, 
yielding after cancelation: 

L[f]<-^{Df,Dv^)g. 

Applying the (weak) maximum principle to /, we conclude: 

max < max f < max f < v"^ max 
E E dj,E' - OpE ' ' 

which implies the result, with an explicit constant oq = vh^. 

Remark. In the general (not necessarily concave) case, we have: 

L[\h\^] < -2|V/i|2 + c„|/i|^ 

Then the proof follows the same lines as that of proposition 3.21 in [2]: we 
apply the maximum principle to / = \h\'^ {"q o v"^) , for a carefully chosen 
function r/(s). 

Evolution o/|V/ip. 

In the calculation that follows, we adopt the usual convention that in 
symbols such as: 

V^/i * (V/i)(2) * h^^\ (V-'7i)(P) = V-'/i * . . . * V^7i {p times ) , 

* denotes some unspecified ^(-contraction of the tensors in question. 
For the time derivative, we have: 

dt\vh\^ = 2{dtivh),vh) + at(5'^/v")(v^/i)pr(v,/i),. 

= 2{dt{Vh),Vh) +3{dtg'^){V,h,V,h), 

using the Codazzi identity. 

For the Hessian (using Vkdi = h^iuj, see appendix): 

V|^J V/i|2 = 2{Vi{VkVh),Vh) + 2{VkVh, ViVh) - hMVh\'' 

= 2{Vli{Vh),Vh) + 2{hkiV^Vh, Vh) + 2(VfcV/i, V/V/i) - hkidjVh 
= 2{Vli{Vh),Vh) + 2{VkVh,ViVh), 
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after cancelation. Taking traces we find: 



{dt-A)\Vhf = -2\V^h\^ + 2{{dt - A){Vh),Vh) +3{dtg'^){Vih,Vjh). 

Commutation of covariant derivatives introduces the Riemann curvature 
tensor, and the time derivative of the connection is also needed: 

{dt - A)(V/i) = V[{dt - A)h] + {VRm) *h + Rm* (V/i) + {dtV) * /j, 

where (see appendix): 

dth = VdH + HV^h + T + hP^ , Tij = Hih{iu, dj) + Hjhiuj, di), 

which combined with T = hcu and dtto = VH + h^"^^ is easily seen to imply: 

dtV = {VdH)u + Vh*h + /i(^) ~ V^^ + Vh*h + h^^\ 

Prom the Gauss equation, Rm ^ h*h. Thus: 

{{dt-A){Vh),Vh) - {V[{dt-A)h],Vh)+V^h*Vh*h+{Vh)^^hh^^hVh*h^^l 
On the other hand, from the evolution equation for h (appendix) we have: 

(v[(9t-A)/i],v/i) = (v(Fv^/i+r+/i(3)),v/i) = (v(i?v^/i),v/i)+(vr,v/i)+(v/i)(2) 

Computing the terms on the right, we find: 

{V{HV^h),Vh) = {V^h,VvHh)+H{V{V^h),Vh) = {V^h,VvHh)+V^h*Vh*h, 
and using the Codazzi identity: 

(VT, V/i) = 2{V^h, VvHh) +V^h*Vh*h + (V/i)^^) * /i^^^ . 
Putting together these results, we have: 

{{dt-A){Vh),Vh) = 3{V^h,V^Hh)+V^h*Vh*h+{Vh)^^Uh^^^ +Vh*h^^\ 

On the other hand, using the expression for dtg^^ given in the appendix we 
find: 

3dtg'^Vih,Vjh) = -6{V^h,VvHh) + {Vh)^^^ * h^^\ 
So we have a cancelation, and obtain the evolution equation: 

{dt - A)|V/i|2 = -2|v2/i|2 + V^h*Vh*h + (V/i)^^) * h^^^ +Vh* h^^\ 
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Remark: without the cancelation, the right-hand side would involve terms 
of type {Vh)^^\ which would be a problem for the argument that follows. 

Given this calculation, the following lemma has a very simple proof. 

Lemma 13.2. For a solution w G C^'^{E), assume we have a uniform 
bound for h: \h\ < oq in E. Then we may find constants a > 0, C > 
depending only on the dimension and oq, so that the function: 

f{x,t) = a\Vh\'' + \h\'' 

is a subsolution in E: {pt — A)/ < C. 

Proof. The calculation above implies: 

A)|V/ip < -2|v2/i|2 + c„(ao|V2/i||V/i| +a^|V/ip + a^|V/i|), 
while the evolution equation for \h\^ implies: 

{dt - A/i)|/i|2 < -2|V/i|2 + Cn{al\Vh\ + 4). 

Clearly we may choose a small enough to satisfy the claim. 

Our next goal is to extend this argument to higher covariant derivatives 
of h. It turns out this does not involve a cancelation similar to the one 
noted above. The various terms T appearing in each expression below are 
all of the same 'weight', where the weight ii; of a term is the positive integer 
defined by: 

^v^7l]=i+l, ^(v^7t)(p)] = + 1), i>o,p>i, 

r 

If r = (V^i/i)(Pi) * ... * (V>/i)(P'-\ w[T] = Y^PiUi + !)• 

1=1 

We introduce a convenient notation for the 'error terms'. For integers wq > 1 
and n >, the notation £''"o,n jg used for a generic term of weight wq, involving 
covariant derivatives of h of order at most n, satisfying certain restrictions: 

w[T] = Wo, ji <n, Pi = lit ji = n, = 1 or 2 if ji = n-1 and n > 1. 

The symbol E'^o,n (denotes the larger space where we drop the restrictions: 

(Sometimes the same notation is used for the real vector space spanned by 
terms of the given type). For example, above we checked that: 

(at-A)|/i|2 = -2|V/i|2 + £;4'\ 
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{dt - A)|V/ip = -2\V^hf + E^'^. 

These symbols have some useful properties. For example, it is easy to prove 
by induction that: 

In the proof of this the following easily verified fact (for n > 1) is used: 

^n+3,n+l ^ (V"+l/i) * /t + (V"/l) * [Vh + + ^"+3>"-l. 

The property verified above for low n holds in general: 
Lemma 13.3. For ra > 0: 

{dt - A)|V''/ip = -2|V"+^/ip + £;2"+4,n+i_ 

Proof (for n > 2.) With the natural multi-index notation: 

dt\Vh\^ = 2{dt{V^h),V^h) + dt{g''gP'gn{^'}hU{V^jh)rs, \I\ = \J\ = n, 

where, using the Codazzi identity and the curvature tensor repeatedly: 

5t(/W)(V?%,(V}/i),, = (n + 2)(5t5'^KV,V"-^/i,V,V'^-^/i) 

+(%'^)i?m[V"-2/i]i * i?m[V"-2/i]j + (%^^)(V"/i)i * i?m[V"-2/i]j-. 

Recall (appendix) dtg^^ = Vh + h^"^^ and Rm = h* h. Thus this term has 
the form: 

(V"/t)(2)*(V/i+/i(2))+(V"-2/i)(^U(V/i+/t(^))*/i(^)+(V"/i)*(V"-^/i)*(V/i+/i(2))*/i(2), 

so it is in £;2"+4,n+i 

Turning to space derivatives, we have (as for n = 1): 

A|V"/i|2 = 2(A(V"/i),V"/i) +2|V"+^/l|^ 

and therefore: 

{dt - A)|V"/ip = -2|V"+^/ip + 2{{dt - A)(V"/i), V"/i) + 
Claim, {dt - A)[V"/i] G E^+^.^^+i (n > 0). 
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The claim clearly implies the conclusion of the lemma (bearing in mind 
the expression seen above for a general term in ^"+3,™+!-) Qf course, we 
prove it by induction on n, the cases n = 0, 1 having already been checked: 

{dt-A)h = HVu,h+T+h^^^ e E^'\ {dt-A){Vh) = V[{dt-A)h]+E'^''^ e E^''^. 

Proof of claim. The claim would follow inductively from: 

{dt - A)[V"+^/i] = V[{dt - A)(V"/i)] + E"+4.«+2^ 
since V£;"+^'"+^ C £;"+4,«+2, 

For the time derivative part, we have, for any multi-index il of length 

i + 1, |/| = n: 

5t[V"+i/i]a = dt[di{v^h[di\)) - (v"/i)(v,9/)] 
= diidtiv^'hidi])) - dt{V'%v^di]) 

= Vi{dt{v^h))[di] + dt{v^h)iVidi) - dtiv^h){Vidi) - v^h[dt{Vidi)]. 

For a multi-index I = ii . . .in of length n, denote by Ip the multi-index of 
length n obtained from I by setting its k^^- entry i^ equal to p. It is then 
clear that: 

k=l P 

In symbolic notation, the preceding calculation is summarized as: 

dtiv'^-^^h] = v(atV"/i) + (V"/t) * (dtV). 

Since dfT e E^'^, this says: 

Covariant derivatives in space may be dealt with in the usual way. Writ- 
ing a multi- index / of length n + 1 as I = il' ,\I'\ = n, we have for first-order 
derivatives: 

Vik(V7+^/i) = VkiVlVh) = Vfe(V,(V?,/i)) - Vfc(V"/t(V,5r)) 

= Vi(Vifc(V7,/t)) - Vik(V"/i(V,5jO) + Rmik[V'},h]. 

And for second-order covariant derivatives: 

VKVfe( V7+^/i)) = VKV,(Vfe(V?,/t))) + Vi{Rmik[V^,h]) - Vi{Vk{V^h{Vidr))) 

= Vi(VKVfc(V7,/i))) + RmalVkiVrh)] + V{Rm * V"/i) + V2(V"/i * h) 

= Vi(V;\(V^,/i)) + Vi(Vv,9,V?,/i) + Rm * V"+^/i + V(i?m * V"/i) + V2(V"/i * h). 



47 



^Iki^r'h) = v,(v2,(v?,^)) - VyM^r'h) + v(r * v-+'h) 

+Rm * V"+i/i + V{Rm * Vh) + V2(V"/i * h) 

= Vi(v2fc(V7//i)) + r * v"+2/j + v(r * v^+^/i) 

+Rm * V"+i/i + V{Rm * V"/i) + V^iVh * h). 
Taking traces with g^^ and using 

Rm = h^'^\ VRm = Vh*h, V = hu, VV = Vh + h^'^\ 

it follows easily that: 

/^{V^+^h) = V(A(V'^/i)) + 

and therefore: 

{dt - A)[V"+^/i] = V[{dt - A)(V"/i)] + E"+4."+2^ 

proving the claim and the lemma. 

The analogue of Lemma 13.2 for higher covariant derivatives of h follows 
easily from these remarks. 

Lemma 13.4. For a solution w € C^^^^'^^~^^'^ {E) , assume we have a 
uniform bound for h and its first n covariant derivatives: \V^h\ < aj in E, 
j = 0, . . . ,n. Then we may find constants a > 0, C > depending only on 
the dimension and the aj, so that the function: 

/„+i(a;,t) = a|V"+^/i|2 + |V"/i|2 

is a subsolution in E: {dt — A)/„+i < C. 

Proof. In the proof we denote by C„ a generic positive constant depend- 
ing only on dimension and the aj, j = 0, . . . ,n. We have: 

(at-A)|V"/ip = -2|V"+^/ip+£;2"+^''*+\ (at-A)|V"+^/ip = -2|V"+^/ip+£;2"+6 
where: 

^2n+4,n+l ^ V"+^/l* V"/l*/l+ (V''/l)(2) + (V'^/l) +^2n+4,n-l ^ 

This implies: 

{dt - A)|V"/i|2 < -2|V"+i/i|2 + C„|V"+^/i| + C„, 
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(5t-A)|V"+i/i|2 < -2|V"+2/j|2+(^^|v«+2/i||v"+i/i|+C7„(|V"+i/i|2+|V"+i/i|+l). 

It is easy to see from these inequalities that we may choose a sufficiently 
small so that the conclusion of the lemma will hold. 

14. Holder gradient estimate for the second fundamental form. 

Notational remark: In this section, parabolic Holder spaces are denoted 
by a single superscript (that is, the space previously denoted (72+q!,(i+q!)/2 
is denoted C^"*"" in this section). Capital X,Y etc. denote general points in 
the spacetime domain E. (This follows the notation used in [8].) 

A continuation criterion for the solution w{y,t) in E'^ in terms of a 
bound on the norm \h\g of the second fundamental form would follow from 
an a-priori C^+\E^) bound on a solution, assuming \ h\g < in E^; equiv- 
alently, from a global a priori Holder gradient bound |V/i|5 < M in E^ (for 
suitably controlled M). In this section we show how such a bound follows 
from the a-priori estimates of linear parabolic theory applied to the evolu- 
tion equations for v, H, and the Weingarten operator, under an additional 
hypothesis. 

Assuming w G C^+^(£''^) is a solution, satisfying in addition \h\g < oq 
m E^, we already observed the maximum principle implies bounds: 

— ~ IT 

< w < Wo, 1 < V < V in E , 

depending only on the initial data and /3 (we assume w > 0, at t = 0, 
vanishing only on ODq.) In particular, g is uniformly equivalent to the 
euclidean metric on E^. In this section, bounds depending on uq, v and the 
initial data will be denoted generically by a constant M > (dependence 
on f3 will not be recorded explicitly). The bound on h implies a uniform 
bound for the spacetime domain E'^ , which we can express in terms of 
a diffeomorphism $ : Dq x [0, T] E'^ by: 

\^c^ < M. 

We will also need to assume a uniform gradient bound on the boundary 
for the second fundamental form: 

\{Vrh){T,T)\<ai Vr G TaZ)(t),|r| = 1. 

Estimates depending oq , ai , v and the initial data will be given in terms of 
constants denoted generically by Mi. 
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In fact E is a bounded domain in W X [0, T] of class (7^+^ with bounds 
controlled by Mi. (This statement includes some regularity in t, so it is not 
immediate from the uniform bound assumed for v*""/!***" on diE). To see 
this, consider the equation satisfied by Wk = d^w, written in 'divergence 
form' with Dirichlet boundary conditions: 

f dtWk - di{g'WjWk) = := {dkg'^)wij - {dig'^)djWk 
I Wk\diE ■= = iPo/P)n'', Wk\t=o = dkWo 

Assuming d^w G C^^^{E), the following estimate holds ([8j, thm 4.27): 

\Wk\l+5 < C{sup\Wk\ + ||/||l,„+l+5 + \ip''\l+5;diE + \dkWo\i+S;Do)- 
E 

Here | l^'^l |i^ri,+i+<5 is the norm in the spacetime Morrey space L^'^^^^^ (E): 
ll/l|i,n+i+5= sup (r-(«+i+^) / \g^\dX). 

Y€E,r<diam{E) J E[Y,r] 

In the present case this can easily be estimated, since: 

\dkg'^\ = \hiJ + h{uj'\ < M, \djWk\ < vao < M ^ |/| < M, 

and \E[Y,r]\ < Cr'^+^ , while 5 E (0,1). Thus | l^^'l < M. 

Since |V^(V^n)| < c(|(V^/i)(t, r)| + \h\) < Mi, it follows that n is in 
space variables on diE. On the other hand, Dw = lo/ (3 on dD{t), and iv is 
a solution of dtu}^ = tVgD'^LO^ + hence n is also in time on diE. 

We conclude \(p%+5;diE < {Po/ P)Hi+5;diE < Mi. 

Therefore we have < Mi, and |?y|2+5 < Mi (note that C de- 

pends on I (7*-' lea and other constants also controlled by M.) In particular, 
E'^ is a C"^^^ domain with chart constants controlled by Mi. (In fact, in a 
neighborhood of any point P G diE with dy^w 7^ 0, a boundary chart ^ is 
given by ^'(yi,y2,t) = {yi,w{y,t),t).) 

The first-order term in the evolution equation for h (or for the Wein- 
garten operator) involves DH; hence the next step is to obtain a global 
gradient bound \DH\i^a 

< Ml in E'^. The mean curvature satisfies the 
'divergence form' equation with Neumann boundary conditions: 

r dtH - d,{g'^{X)Hi) + dj{g'^){X)diH - c{X)H = 0, 
\ dnH = {(3y(3o)HKn := ^ on dD{t), H\t^^ = Hq, 

where: 

c := \h\l - h^{LO,Lu) + Hh{uj,uj). 
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Then, with the regularity conditions for the domain and the coefficients: 
and assuming H G C^'^^{E), or w e C^'^^{E), we have the bound: 

\H\i^S.E < C{sup\H\ + \ip\s;diE + \Ho\i+S;Do)- 
E 

As noted ear her: 

\\djg'^\\l,n+l+5 + ||c||i,„+i+5 < M, 

hence C is controlled by M. In addition, |io|2+5 < Mi implies \h\g < Mi, 
and hence: \'^p\s■^^E < Mi. We conclude < Mi, and state it as a 

lemma. 

Lemma 14.1. Let w G C^^^{E'^) be a classical solution of graph mean 
curvature motion, with contact and constant-angle boundary conditions. 
Assume (i) \h\g < oq on diE; (ii) |(Vi-/i)(t, t)| < ai on diE. Then we have 
the global gradient bound for H: 

sup \DH\s < Ml, 

ET 

for a constant Mi depending on 5, v, cq, ai and the initial data wq. 

Corollary 14.2. Under the same hypotheses as Lemma 14.1, we have 
a global gradient bound: 

sup|V/t|p < Ml, 

E 

for a positive constant Mi depending on S,v,ao,ai and the initial data wq. 

Proof. The bound on the components (\/nh){T,T) and (SJ nh){n,n) on 
the lateral boundary diE follows immediately from the expressions in section 
11. The bound on (VT/i)(T, r) over diE is hypothesized, and then the bound 
on the remaining component {'Vrh){n,n) follows from the global gradient 
bound \DH\ < M implied by Lemma 14.1. Thus |V/i| < Mi on diE, and 
then the global bound follows from Lemma 13.2 and the maximum principle. 

To improve the conclusion of corollary 14.2 to a Holder gradient bound, 
it is natural to consider the evolution equation for h, with the Neumann- 
type boundary conditions derived in section 11. One is then faced with the 
problem that those boundary conditions do not control components such as 
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(Vr^)(T, r) on dtE. So as a preliminary step we consider the evolution equa- 
tion for f , which has the advantage that the boundary values are constant. 
Written in linear form, we have: 

dtv - g'^{X)vij + b\X)diV + c{X)v = 

V\dtE = ^, V\t=0 = VQ, 



where: 



^^'^''^ = ^^^-iTW^^^' ^^^^^ = ilW^'')' ^(^) = l^l^(^)- 

We clearly have g^^ G (since Dw G C"'), as well as 6*,c G (since 
h e C^) and diE e C'^'^^, with bounds controlled by Mi in all cases as 
observed earlier. Therefore assuming v G C'^+^ (equivalents, w G C^+^) we 
have the bound: ^ 

\v\2+S;E < C(supi; + -), 
E P 

with C controlled by Mi. Thus |-D^?^|5 < Mi. Recalling v~^diV = h{di,uj), 
this implies: 

\{Vrh){n,n)\5-dtE = \{Vnh){T,n)\5-AE < Mi, Vr G T9D(t),|r| = 1. 

Since H = (S'^hnn + h{T, r) on diE, it follows from Lemma 14.1 that we also 
have \(Vrh){T,T)\s-djE < Mi. For the remaining components of V/i, this 
bound follows directly from the boundary conditions: 

\{Vnh){T,T)\s;d^E + \{^nh){n,n)\s;a,E < Ml. 



Now consider the evolution of the components of the Weingarten oper- 
ator, written in divergence form with Neumann boundary conditions: 

J dth'^ - diifdih]) = in E'^, := HjhfJ - Hih]u^ + hf^ - {di(f^){dih]) 
\ dn{h]) = ^) on diE, h)^^^^ = h% 

The same theorem quoted above gives the estimate (assuming G C^^^ , 
or we C^+^): 

\h%+5;E < C{sup\h';\ + ||/j=||i,„+i+5 + \^';\S-AE + \h'^j\ l+S;Do)- 
E 
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Note that: 

dnih';) = g'^{Vnh)ij = P^{Vnh){n,dj)n'' + {Vnh){T,djy on diE. 

Prom this and the above discussion it fohows that \y^j\s;diE ^ Mi. The 
bound I l/j^l |i,n+i+<5 ^ Ml follows from Lemma 14.1 and Corollary 14.2. We 
conclude \hJ\l_^_§.^ ^ -^i- The 1 + 6 estimate for hj clearly implies the 
following lemma. 

Lemma 14.3. Let w G C^'^^{E'^) be a classical solution of graph mean 
curvature motion, with contact and constant-angle boundary conditions. 

Assume (i) \h\g < ciq on diE[{n)\{Vrh){T,T)\g^E < cli- Then we have the 
global Holder gradient bound for h: 

\Vh\s.ET < Ml, 

for a constant Mi depending on 5, v, qq, ai and the initial data wq. 

Remark. This is clearly equivalent to a global 'a priori' C^'^^ bound for 
w on E'^, \w\s-^-s < Ml. 

Lemma 14.3 is the main step in the derivation of a 'continuation criterion' 
for this flow. 

Proposition 14.4. Assume the maximal existence time T^nax is finite. 
Then (for n = 2, in the concave case): 

limsup sup {\h\g + |(Vr/i)(T, t)|) = oo. 

t^Trr^ax dD{t) 

Proof. For wq G C^'^°'{Do) satisfying the contact angle condition (with 
a e (0,1) arbitrary) and a = 6?, Theorem 8.1 yields a unique solution 
F = [u, v?] of m.c.m. with contact angle/orthogonality boundary conditions 
in a maximal time interval [0,T^aa;], with F G C2+"((5o'"""), Qg""^^ = 
Q X [0,r^aa;); this is also the unique solution in F G C^+^^ (Qq'"'"^), where 
(5 = o? . Then w = uo (p~^ G C^+^(£^'^'""^) is a solution of graph m.c.m., 
which for any to > is in C^~^^ (eJ^'"'^). By contradiction, assume + 
I {'Vrh){T, t)| is bounded in Ef^ for any T < Tmax (with bound independent of 
T). Then lemma 14.3 applies, giving an a-priori bound |V/i|r.^T < Mi, for 

T arbitrarily close to Tmax- In particular, \w{-,T)\^a+s(^ij^j'y^ < Mi, and for 
T sufficiently close to Tmax we can use theorem 8.1 again (with initial data 
w{-,T)) to find a solution F' = [u',^'] G C^+^'iQ^') (where T' > Tmax), 
extending F. This contradicts the maximality of Tmax- 
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15. Behavior at the extinction time. 

In this section we consider the behavior of Sf as t approaches the maxi- 
mal existence time T, in the concave case. We also assume H < Hq < at 
t = 0, so T is finite. Let Kt C R""*"^ be the compact convex set bounded by 
Sf. Since H < 0, {Kt} is a decreasing family, and the intersection: 

Kt= Pi KtC M"+^ 

0<t<T 

is compact, convex and non-empty. It turns out Kt has zero (n-|- l)-volume. 
In this section we use the support function to show this when n = 2 (follow- 
ing the argument in [l3]), under the assumption that there is no 'gradient 
breakdown'. 

Assume the origin S R"'^^ is a point of Kt- The support function of 
Kt (with respect to this origin) is the function p{-,t) on D{t): 

p{y, t) = {G{y, t),N{y, t)), G{y, t) = [y, w{y, t)]. 

Since Kt is convex, p > in D{t); the evolution equations and boundary 
conditions for p are easily computed: from L[G] = and L[N] = \h\gN, we 
have: 

L[p] = {L[G],N) + {G,L[N]) - 2g^\dkG,diN) 
= \h\lp + 2H, 

and since {dnG, N) = 0: 

PnldDit) = {G,dnN) = -A{G^ ,N), 

where, with y'^ '■= y — {y ■ n)n G TydD{t), the tangential component G^ : = 
G — {G, N)N is easily seen to be, at dD{t): 

G^ = lKy^ + y,0]. 

Since h{y'^,n) = at dD{t), this implies A{G'^,n) = j3'^{y ■ n)h{n,n). Note 
that p{y) = —(3o{y ■ n) on dD{t), so we have: 

Pn\dD(t) — -rPiT-nn, 

m 

reminiscent of the boundary condition for H. Note also that we have the 
upper bound: 

p<\\G\\ < max||Go|| := Po, 
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since the Kt are 'nested'. 

Proposition 15.1 (n = 2) Assume we know that 

Umsup( sup \h\g) = CO 

t-^T y€dD(t) 

at the maximal existence time T. Then: 



Hminff inf p(y,t))=0. 



Proof. Reasoning by contradiction, assume p > 2(5 > for i G [0, T) . 
We claim that this implies an upper bound for \H\ (and hence for |/t|, since 
< nH"^), contradicting the fact that lim sup^_,.2- supp^ \h\ = oo. 

To prove the claim, consider the function: 



p — 5 p — S 

Using the evolution equations and boundary conditions for H and p we find 
(with a) := Lo/\uj\g, see Remark 15.1 below): 

m = fi-ml + 2p/ + \u;\l{h\co,co) - Hh{u,^)) - -^^g'^^dkfdip, 



p — 



/n|aD(t) - -^J^hnn-^^ > 0. 



Since h\l > (l/n)i72 = {l/n)f{p-Sf: 

L[f] < JP-^^'^ f + 2pf+\u\l{h\u;,u;) - Hh{u,u)) - -^{Vf,Vp)g. 

Now recall (Remark 12.1) that if n = 2: 

h^{u;,6j) - Hh{u,u) = -A < 0, 

so: 

Assume (5 > is so small that sup^^Q-j /|t=o < ^f?^- We claim this persists 
for all t G [0, r). If not, assume f{yo,to) = 2npo/5^, with to > smallest 
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possible and yo a local max of /(-jto)- Since /„, > at dD{tQ), zq = (yoj^o) 
can't be a boundary point of E (by the boundary point lemma). Since 
yo G dDito) is an interior point, we have -^^[/]|zo ^ ^fi^o) = 0) hence: 

-f[zo) < 2p[zo), or f{zo) < y? Jy? < 



which is not possible since p — 6 > 6. Thus f{y,t) < Apo/S^ in which 
implies the bound: 

for t G [0, T), contradicting the maximality of T. 

Remark 15.1. Consider the vector fields in D{t) C M^: 

UJ = -[wi,W2], = VUJ'^ = \ — W2,Wl]. 

V 

It is easy to verify the following: 

{uj,uj)g = 0, |a;|^ = ItDl^ = \Dw\l := wl + wj. 

Thus we may think of {to, to} as a 'conformal pseudo-frame' (w and oo vanish 
when Dw = 0), defined on all of D{t). Moreover, at the boundary dD{t): 

w = ii^n, UJ = —n := — -r, 

f3 (i 

where {t, n} is an euclidean-orthonormal frame along T^. Thus to and Co 
supply 'canonical' extensions of n, r to the interior of -D(t), as uniformly 
bounded vector fields. 

Remark 15.2. It follows from the proposition that Kt cannot con- 
tain a half-ball of positive radius centered at a point of M^; in particular, 
vo13{Kt) = 0. Based on the experience with curve networks [7J, one is led to 
expect that Kt is a point (that is, that diam{KT) = 0), at least under the 
same assumption as the proposition ('no gradient breakdown'). We have not 
been able to show this yet; existence of self-similar solutions and comparison 
arguments appropriate to the free-boundary setting appear to be needed for 
the usual approach to work. 

16. Final comments. 
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1. We state here the local existence theorem for configurations of graphs 
over domains with moving boundaries. In this setting, a triple junction 
configuration consists of three embedded hypersurfaces in M"^"^, 
graphs of functions defined over time-dependent domains D^(t),D'^{t) C 
M" [D^ covered by one graph, by two graphs), satisfying the following 
conditions: (1) The T,^ intersect along an (n — l)-dimensional graph A{t) 
(the 'junction'), along which the upward unit normals satisfy the relation: 
-^1 + -^2 = -^3- (2) If a fixed support hypersurface S C is given (also 
a graph, not necessarily connected), the intersect S orthogonally. 

Topologically, in the case of bounded domains one has the following 
examples: (i) ('lens' type) 2 disks (or two annuli) covering D'^{t) and one 
annulus covering D^{t); (ii) ('exterior' type) two annuli covering -D^(t) and 
one disk covering D^{t). The boundary component of the annuli disjoint 
from the junction intersects the support hypersurface S orthogonally for 
each t. 

Let El (J = 1,2,3) be graphs of C3+" funct ions over C^~^" domains 
Dq,Dq C M", defining a triple junction configuration and satisfying the 
compatibility condition for the mean curvatures on the common boundary 
To of and D^: 

+ = H^. 

Then there exists T > depending only on the initial data, and functions 
g (^2+a,i+a/2(Q/)^ c M" X [0,r), SO that the graphs of w^{.,t) : 

D\t) R define a triple junction configuration for each t S [0,T), moving 

by mean curvature. 

The proof will be given elsewhere. 

2. An interesting issue we have not addressed here is whether one has 
breakdown of uniqueness for initial data of lower regularity, or if the 'or- 
thogonality condition' at the junction is removed. For curve networks, non- 
uniqueness has been considered in [llj ; but neither a drop in regularity (from 
initial data to solution, in Holder spaces) nor the orthogonality condition 
play a role in the case of curves. 

Appendix 1: Proof of lemma 4.1. 

Throughout the proof, n denotes the inner unit normal at dD, extended 
to a tubular neighborhood M so that D^n = 0. Since D is uniformly C^"*"", if 
follows that n S C^^"(9D), with uniform bounds. Denote by p the distance 
to the boundary (so Dp = n in M). Let C ^ C^{D) be a cutoff function, 
with C = 1 in M C AA, C = in D \ AA. 
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We find ip of the form: 

(^(x) = X + ({x)f{x)n{x) 

with / G C'^^°'{M). The 1-jet conditions on (p at dD translate to the 
conditions on /: 

f\dD = 0, Df\Qo = 0, D^f{n, n)i9D = ^f\QD = h. 

Now use: 

Lemma A.l. Let I? be a uniformly C^+" domain with boundary distance 
function p > 0. Let h G C'^{dD) be a bounded function. Then there exists 
an extension g G C°°{D) PI C{D) so that g^g^ = h, sup^ \g\ < supg^, \h\ and 

Given this lemma, all we have to do is set / = (l/2)p'^g, which clearly 
satisfies all the requirements (in particular, A/ = h at dD.) 

To verify that is a diffeomorphism, it suffices to check that \Cfn\ci (in 
J\f C {p < po}) is small if po is small. This is easily seen: 

\Cfn\co < {l/2)pl\g\co; 
\DC\<cpo^ ^\fDC\<cpo\g\co. 

\Df\<il/2)pM C2+«(D) 

on A^, since Df G C^+"(L>) and Df^gu = 0. And finally, with A the second 
fundamental form of dD: 

\Dn\ < \A\co ^ \fDn\ < il/2)pl\g\co\A\co. 

A word about Lemma A.l. (This is probably in the literature, but I 
don't know a reference.) If D is the upper half-space, we solve = 0'm D 
with boundary values h. Then the estimate 

[D\p''P*h)]^°'\D)<c\h\c^^SD) 

follows by direct computation with the Poisson kernel P; for the rest of the 
norm, use interpolation. Then transfer the estimate to a general domain 
using 'adapted local charts', in which p in D corresponds to the vertical 
coordinate in the upper half-space. (It is easy to see that at each boundary 
point there is a C^"*"" adapted chart, with uniform bounds.) 
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Appendix 2: Evolution equations for the second fundamental form. 

We consider mean curvature motion of graphs: 

G{y,t) = [y,w{yM yGl?(t)cM", 



wt = g^-^Wij = vH, V = ^1 + \Dw\'^. 

In this appendix we include evolution equations for geometric quantities, in 
terms of the operators: 

dt -Ag, L = dt- trgcf. 
It is often convenient to use the vector field in D{t): 

LO := -Dw. 

V 

Since —u> is the component of the unit normal N and L[N] = \h\^N, we 
have: 

L[u'] = \h\y, \h\l:=g'^g^%,hki. 
Here h = (hij) is the pullback to D(t) of the second fundamental form A: 

h{di,dj) = hij = A{Gi,Gj) = ^Wij. 

First, denoting by V the pullback to D[t) of the induced connection 
(that is, G^iVxY) = '^g,xG*Y for any vector fields X,Y in D{t)), and 
using the definition: 

V^.Gj = Gij~{Gij,N)N = [Q,Wij]-^Wij[-Dw,l] = ^[Dw,\Dw\^] = ^G,Dw, 

we conclude: ^ 

V didj = —hijDw = hijio. 



Prom this one derives easily a useful expression relating the Laplace- 
's' 



Beltrami operator and the operator trgCp acting on functions: 



^gf = trgcff - —Wmfm = trgd^f - Hdu,f. 



We also have, for the covariant derivatives of h with respect to the eu- 
clidean connection and to V = V^: 
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(Here is the symmetric (3, 0)-tensor with components: Vmhij = (Va^/t)(Oi, Oj).) 

Iterating this and taking ^-traces yields (using the Codazzi identity and 
the easily verified relation diUj'^ = := g^'^hij): 

tr,d2(/i,,) =5™'=a„(afc(/i,,)) =/"'=(Vi„^,,,^/i)(5,,9,) 

+HV^hij + 2[h^Vkhjp + h^Vkhip]ujP + [Hihjp + Hjhip]^^ 
+2[hip{h^)jq + ih%phjq + HhiphjqjujPuj'i + 2{h%j + 2{h%jh{uj,cj). 

Here the powers hP' and of h are the symmetric 2-tensors defined used 
the metric: 

{h% := g^^hikhpj = h^hpj, {h% := g''^ g^" hikhpihqj . 
Note also that: 

[hi^khjp + h]Vkhip]uP = Vu,{h%, 
using the Codazzi identity. 
Evolution equations for h. 

Starting from Gt = vHe-a+i = H{N + l[Dw, \Dw\^]) = HN + HG^uj 
and Nt = -V^H - Hv^V^v (where V^/ = g'^ fjGi and V/ = g'^ fjdi) we 
have: 

dt{hij) = {{HN)ij,N) - {Gij,V^H) - ^{Gij.V^v) + {{HG.u)ij,N). 
Using the easily derived facts: 

{Ni,,N) = -h^{^^,^j), 

Hij - {Gij,V^H) = {VdH){di,dj), 
^{Gij,'V^v) = h{uj,uj)hij, 

we obtain: 

dtihij) = {VdH){di,dj) - Hh^{di,dj) - Hh{co,iu)hij + {{HG,u>)ij,N), 
where: 

{{HG,Lo)ij,N) = Hi{{G,uj)j,N) + Hj{{G,u)i,N) + H{{G,u;)ij,N). 
To identify the terms, computation shows that: 

{{G,u;)i,N) = h{LO,di), 
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and hence, using also: 

V§.(G*u;) = G*(Va.u;), Vg.o; = (/if + u;'^hiquP)dp = ^ hipdp, 

p 

we obtain (using uj^dj{hik) = ^uhij + 2h{di,uj)h{dj,uo)): 

{{G^uj)ij,N) = dj{uj''hik)-{V^.iG^u;),djN) = h]hik+u^dj{hik)+h{dj,V Si^) 

= iSujh)ij + {h?)ij + 2h{ijj, di)h{uj, 3j) + ^ hiphjp 

p 

= {^coh)ij + 2{h%j + 3/i(a;, di)h{u, dj), 

since J^p^ip^jp — + h{oo,di)h{uj,dj). Combining all the terms yields 

the result: 

dtihij) = {VdH){di,dj) + HV^hij + Hih{u;, dj) + F,/i(w, di) 
+H{h^)ij + SHh{u}, di)h{u>, dj) — Hh{u>, u>)hij. 

Prom this expression and Simons' identity (in tensorial form): 
VdH = /:^gh + \h\lh- Hh^, 
we obtain easily a tensorial 'heat equation' for h: 

[{dt - ^g)h]ij = HVojhij + Hih{uj, dj) + Hjh{u, di) + Cij, 
Cij := \h\ghij + 3Hh{di,u)h{dj,u>) — Hh{uj,uj)hij. 

Using the earlier computation relating Agh (the tensorial Laplacian of 
h) and trgCph, we obtain from this the evolution equation in terms of L: 

L[hij\ = -2[h'tV^hjk + h]V^hik] + Cij, 

Cij ■= Cij-2[h{di,u)h'^{dj,u)W{di,u)h{dj,u)]-2{h%j-2{h%jh{u,u)-2Hh{d^ 
We may also write this purely in terms of the euclidean connection d: 

L[hij] = -2[hfdu;hjk + h^d^hik] + Cij, 

Cij = Cij+2[h{di,u)h^{dj,i^)W{di,u)h{dj,u)]-2{h%j-2{h^)ijh{u,u)-^ 
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Time derivatives and evolution equations for uj and g. 

It is sometimes convenient to use the 'Weingarten operator': 

S{X) := S{X'di) = h)X^di. 

The time derivative of u is simply minus the time derivative of the 
component of N. In addition, one computes easily that ^ = S{ijo), so we 
have: 

dtuj = VH+ —Vv = VH + HS(lo). 

V 

For the metric and 'inverse metric' tensors we have: from dtgij = {wiWj)t 
and Wit = {vH)i: 

dtgij = v'^{HiJ + HjJ) + v'^H{h{u, di)u^ + h{u, dj)J), 

and then, using dtg''^ = -g'^'^dtgkig''^ ■ 

dtg'^ = -[{VHYu^ + {VW)J] - H[S{uyuj^ + S(uyJ]. 

Since we know the evolution equation of it is easy to obtain that of 

gii: 

L[g^i] = -L[lo'lo^] = -L[u']uj^ + 2g^\dkJ){diu^) - JL[u^. 
Using dkUj'^ = Kf., we find: 

L[g'^] = -2\h\lJuj^ + 2{h'^yK 

It is also easy to see that d^g^^ = —{h],ijj^ + hif^u'^). 
Evolution of mean curvature. 

To compute the evolution equation for H = g^^hij, we just need to 
remember g'^^ is timc-dcpcndcnt: 

idt-Ag)H = {dtg'^)ihij)+trg[{dt-Ag)h] = -2h{V H,u)-2Hh\u,u)+trg[{dt-Ag)h]. 
The result is: 

{dt - Ag)H = Hd^H + \h\lH + Hh^{io,io) - H^h{io,io). 

Since L[f] = {dt — A.g)f — Hd^^f (for any /), we see that the equation in 
terms of L has no first-order terms: 

L[H] = \h\lH + Hh^{uj,uj)-H^h{uj,uj) 
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Remark. One can also find L[H] starting from the expression: 

L[g'^h,] = L[g^^]h,,+g^^L[K,] - 2g^\dkg''){dih,,). 
This may be used to check the calculation. 
Evolution of the Weingarten operator. 

The tensorial Laplacian of S is the (1,1) tensor A^^ with components 
Aghj. We have: 

Agh^ = g^'Agh,,, or {{AgS)X,Y)g = {Agh){X,Y). 

The evolution equation is easily obtained: 

{dt - Ag)h] = {dtg''')hij + g'\dt - Ag)hij 

= HV^h]+Hjh^u^-Hih]u^^h\]h]+2HS{ufh{u, dj)-Hh{u,u)h)-Hh{S{u), dj)u^. 
Remark: Since the components of VS are given by: 

we see that upon setting j = k and adding over k we recover the evolution 
equation for H. 

The evolution equation for h'j in terms of L follows from the calculation: 

L[h'^] = L[g''']h,,+g^'^L[h,,] - 2g^^ {dm9'''){dnh,) 
= -2{V^hl)h^ + {d,\h\l)u:^ 
+\h\]h^^ - Hh{u,u)h'; + HS{Lo)^hidj,io) + 2h^{dj,u;)io'' - 2{h'^)^ujP h{dj , u) . 
Setting j = k and adding over k, we recover the earlier expression for L[H]. 

Evolution of \h\g. 

The fact that g'^^ is time-dependent introduces an additional term in the 
usual expression: 

{dt - Ag)\h\l = -2\Vh\l + 2{h, (dt - Ag)h)g + 2{dtg'^ ) {h^j . 

Using the expressions given earlier, one easily finds: 

{dt - Ag)\h\] = -2\Vh\] + Hd^\h\l + 2\h\^g - 4Hh^{LO,io) - 2H\h\lh{uj,u), 
L[\h?g] = -2|V/i|2 + 2\htg - AHh\u,uj) - 2H\h\lh{u;,u). 
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